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ABSTRACT 


An analytic study was carried out to determine the 
applicability of the concept of thermal stress fragmentation 
to the UO./sodium Fuel-Coolant Interaction. Major emphasis 
was put on the fracture mechanics approach to assessing 
whether or not the solidifying UO, would fracture under the 
thermally induced stresses. It was found that the stress 
levels were sufficient to generate K. values substantially 
in excess of the UO, fracture toughnéss, K Cc: Thus, rapid 
instantaneous propagation of inherent flaws°is anticipated. 
Parametric studies in which the surface heat transfer 
boundary condition was varied did not alter this conclusion. 


Extension of the thermal stress fracture concept to Al.O 
resulted in Similar behavior. Subsequently, this material 
was selected as a good simulant for use in experimental 
studies of the mechanism. Additional studies on tin led to 
the conclusion that thermal stress fracture was not an 
applicable mode of fragmentation inducement in ductile 
metal/water interactions. 


It was concluded that thermal stress fragmentation is a 
feasible mode of fragmentation in the UO.,/sodium system. 
However, further experimental and analytical work is 

necessary to establish firm verification of this model. 
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CHAPTER ONE 


INTRODUCTION 


Possible accident situations, such as the melt down 
of a reactor core with the subsequent interaction of the 
fuel with the coolant, have been of major concern since the 
inception of reactor systems. It is known that sudden 
contacting of a hot liquid with a cold vaporizable one may 
result in a violent explosion. This process, regardless 
of whether explosive or not, will herin be referred to as 
a fuel-coolant interaction (FCI). The hot liguid is the 
fuel, the cold liquid the coolant, and the excess heat in 
the fuel provides the energy which drives the interaction. 
This type of interaction process is not limited to reactor 
systems. However, the fact that toxic radioactivity is 
coupled to the problem in reactors increases the hazards 
involved. Therefore, extensive safety analyses are a 
necessity in an effort to develop design constraints capable 
of minimizing the probability of such an accident occurring. 
These safety analyses will also provide an advance 
assessment of post-accident system behavior, and thereby 
allow for pre-planned measures of damage containment. This 
present work is linked to a portion of this type of 


eZ 
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safety analysis, specifically concerning the liquid metal 
fast breeder reactor (LMFBR). In this particular case, 
the fuel is uo, and the coolant is liquid sodium. 

The specific type of accident referred to above results 
from some type of power excursion, causing fuel element 
failure, which allows the dispersion of molten fuel in the 
liquid sodium. This accident situation, although highly 
unlikely, must be considered. The main concern is whether 
or not such an interaction of large quantities of molten 
fuel and coolant will yield vapor production fast enough to 
result in rapid pressurization coupled with the development 
of a local shock wave. 

Although the precise outcome of the FCI is subject to 
speculation at present, existing models have resulted in 
predictions of sizeable amounts of mechanical work being 
generated from the thermal energy input (49, 65). Normally, 
in these models, a relatively large amount of fuel is 
assumed to interact with a small, constrained volume of 
coolant. Results substantiating such predictions have not 
been obtained experimentally for molten uo. in sodium. In 
one instance, though, where small amounts of liquid sodium 
were injected into molten UO, (the inverse experiment) 
high energy conversion rates were observed (20). Thus, a 
better assessment of the contact modes, energy transfer 
rates between the hot molten fuel and the coolant are 


required to better understand the interaction. 
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For explosive vapor generation to result from a fuel- 
coolant interaction, either an extraordinary heat transfer 
rate must occur, or the contact area for heat transfer must 
be increased substantially. The precise sequence of events 
and the controlling factors leading to such an explosion 
are subject to some controversy. There are essentially two 
basic themes as to what takes place. One centers around the 
spontaneous nucleation theory (74). In this case, perfect 
liguid-liquid contact, a lack of nucleation sites allowing 
for coolant superheating, and the elevation of the interface 
temperature to the spontaneous nucleation temperature are 
required. Spontaneous nucleation is predicted to result in 
extremely high pressure pulses. This mechanism does not 
occur in the case of the UO, in sodium interaction since the 
interface temperature remains below the spontaneous 
nucleation temperature of sodium. 

The other vroposed method of initiating explosive vapor 
generation is based on the assumption that one of the fluids 
(normally the fuel) fragments, rapidly increasing the surface 
area available for heat transfer. Experimental fragmentation 
of uO, in sodium has resulted only in low order 
pressurizations (20). However, if the interaction is 
constrained to a small volume of the coolant as mentioned 
above, existing models predict that potentially large 


amounts of mechanical work could be generated. 
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It 1S not known whether fragmentation is a required 
condition for vapor explosions. Yet knowledge of the 
process itself would permit evaluation of the known 
low-pressure interactions, and the resultant impact on the 
reactor system as a whole in an accident situation. 

The present lack of understanding of the fragmentation 
process ina UO.,-sodium FCI and the inability to accurately 
predict the subsequently generated heat transfer area 
prohibits formulation of an accurately descriptive accident 
model, based on the second theme, for use in design. As 
long as this lack of understanding persists, the time 
dependent, mixing phenomena will remain unknown. 

The broad spectrum of experimental and analytical work 
reviewed in the next chapter demonstrates the wide 
variety of concepts considered thus far. It appears that 
for the specific case of fragmentation in the UO,-sodium 
FCI, the solid shell theories are the most promising. 
Specifically, the thermal stress mode of fragmentation 
during solidification, which is the one under consideration 
here, 1s the most promising. 

This work is primarily intended to further evaluate 
and extend the thermal stress fragmentation theory. Of 
orimary interest is the determination of the solid uo, 
Shell's susceptibility to fracture under the developed 
stresses. Of secondary importance is an evaluation of the 


effect on the fracture potential of gases covering the 
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surface of the uO. particles undergoing thermal stress 
fracture. As well, an estimate of the constraint provided 
by the solidifying UO. surrounding liquid sodium is made 
(the inverse problem). Finally, the applicability of 
thermal stress fracture to other simulant systems is 
included. 

In conducting this analysis, extensive use was made of 
the existing experimental works and analytical models 
described herein. Hopefully, the results advance the 
"state-of-the-art" one step closer to allowing for a 
realistic evaluation of potential LMFBR accident situations 


involving a fuel-coolant interaction. 





CHAPTER TWO 


REVIEW OF EXPERIMENTAL AND THEORETICAL STUDIES OF 


MOLTEN MATERIAL/COOLANT THERMAL INTERACTIONS 


2.1 Introduction 

The complexity of the fuel-coolant interaction and the 
lack of understanding associated therewith are illustrated 
by the extenSive number of published articles relating to 
this and its associated problems. The following is an 
attempt to enumerate and briefly summarize some of the major 
experimental studies conducted to date. Reactor system 
safety analyses, as well as work in non-nuclear fields, 
have generated numerous studies of thermal interactions 
involving a wide variety of materials. Until recently, 
published work concerning molten UO. and sodium was 
somewhat limited. In today's energy-conscious world, 
however, the push for LMFBR completion has stimulated 
a more concentrated effort in this particular area. This 
review is by no means all-inclusive, but the intent is to 


point up some of the important results recorded thus far. 


2.2 Experimental studies 


Active experimental programs have been undertaken 


internationally involving numerous materials. In recent 
Ne) 
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years, reported experimental work related to the molten 
fuel-sodium interactions has come out of the international 
CREST meetings at Grenoble in 1972 and Ispra in 1973, in 
particular. A summary of presentations at the former can 
be found in Ref. 1 and of those at the latter in Refs. 2 
and 3. 

Earlier work concerning the fuel-coolant interaction, 
i.e., GAuring the 1950's and 1960's, was primarily directed 
at thermal reactors utilizing a water coolant. The emergent 
problem was found to be extremely complex and the actual 
mechanisms involved, as well as the interdependence of such 
mechanisms still remain somewhat in question. The majority 
of the work done with respect to the thermal reactor systems 
is not directly applicable to the LMFBR. However, it has 
been the valuable experience gained from such studies that 


has paved the way for the more recent fast reactor analyses. 


2.2.1 Accidental thermal explosions 

Accidental thermal explosions encountered in industry 
resulting from the sudden contact of molten metals with 
damp surfaces or cool liquids have been recorded (4, 5) and 
reviewed (6) covering a time frame extending from the middle 
1800's to the last decade. Some of these accidents have 
resulted in extremely heavy damage, and in some cases have 


even been responsible for numerous deaths. 
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An early indication that a metal-water thermal 
interaction 1S a potential safety problem in reactor systems 
resulted from the core meltdown accident of the NRX Canadian 
test reactor in 1952. This accident, as well as others 
which can be explained in terms of high energy thermal 
interactions are directly responsible for the research push 
into this particular area. Reference 7 provides a good 
summary of these accidents and their causes. Following the 
NRX accident, a set of destructive reactor experiments was 
conducted in the Borax reactors. As reported by Dietrich 
(8), the melting of the fuel plates and the following 
contact of the molten metal with the water coolant resulted 
in the generation of extremely high pressures of approximately 
6-10,000 psi. 

In another accident involving an experimental reactor, 
the EBR-1, the first case of a core meltdown with a coolant 
other than water was encountered. In this instance, the 
coolant was NaK and a large destructive explosion did not 
take place. During post-accident analysis the remains were 
found to be in a porous condition. It was subsequently 
postulated by Kittel (9) that vaporization of entrained 
NaK may have been responsible for the observed spongy 
remains. Post-accident analysis of the well-known SL-1l 
boiling water reactor accident (the first fatal reactor 
accident) pointed strongly to a violent molten metal-water 


reaction following fuel element meltdown as being responsible 
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for the extensive damage incurred (10). The magnitude of the 
pressures developed wasS enormous as the reactor vessel rose 
approximately nine feet off its foundation. Further 
destructive experimental test results, such as those from 

the SPERT reactor tests, tend to support conclusions that 
thermal interactions, and not nuclear explosions, resulting 
from extremely rapid vapor generation cause the destructive 


mechanical work. 


2.2.2 Laboratory experiments 

As a result of the accidents enumerated above, 
extensive experimental research into the mechanisms 
responsible for thermal interactions wasS initiated. Some 
of the earliest experiments were conducted by Long (6), 
followed later by Elgert and Brown (ll). These experiments 
all centered around Molten metal-water interactions. A 
Significant outcome of such studies was the ruling out of 
chemical reactions as a causative mechanism. This hypothesis 
was later verified analytically through a study of chemical 
reaction kinetics by Epstein (12). 

The experiments conducted by Long (6), which are 
becoming something of a clasic in the field, consisted 
of pouring fifty pounds of molten aluminum into a tank 
containing water. He concluded that a thin layer of water 
trapped between the molten metal and the container vaporized 


and thereby initiated fragmentation of the molten aluminun. 
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The resulting fragments allowed an increased rate of heat 
transfer to the surrounding water and thereby created an 
explosion. This mechanism has since been entitled the 
entrapment theory. In an attempt to further substantiate 
his theory, Long demonstrated that when the container was 
coated with grease, no longer allowing entrapment, the 
interaction was suppressed. Further verification of Long's 
results was subsequently provided by Hess and Brondyke (13). 
Other experiments, such as those of Sallack (26) on the 
quenching of paper smelt, culminated in a proposal that the 
mechanism for spontaneous fragmentation of the molten 
material resulted from the vaporization of liquid entrained 


in cracks in the melt's solidified surface. 


2.2.3 Small-scale dropping experiments 

Experimental investigation of the fragmentation of 
small molten drops in cold liquids has been a fairly 
widespread method of analysis. More emphasis has been 
placed on this particular mode of contact as it is believed 
to be somewhat representative of the type of results 
expected ina UO.-sodium free contact situation. A good 
review of the results of these particular experiments can 
be found in Ref. 14. 

Experimenters at Argonne National Laboratory (ANL) 


have conducted numerous dropping experiments consisting of 


dropping various materials into liquid sodium (15-18). Some 
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of the earliest of these experiments were carried out by 
Swift and Pavlick (16). A significant observation was the 
fact that initially solid materials, when dropped, were 
simply quenched by the sodium. Molten samples, on the other 
hand, fragmented, and it was later found that this applies 
to molten uranium (17) and uo. (20), as well. Other 
important results include the decreasing size of fragments 
with increasing initial temperature (18), and the enhancement 
of fragmentation with increased sodium pool temperature up 
to a peak value, followed by decreasing fragmentation with 
further increases (15). 

Resulting from a strong desire for "observability", 
water has become widely utilized as the coolant for 
experimental applications. Swift and Pavlik (16, 17) 
extended the scope of their experiments by using water as 
well as sodium in the coolant pool. A major finding was 
the fact that while some molten metals will fragment in 
sodium, they will not in water. To ensure that these 
results were not dependent on an oxide film effect, they 
utilized noble metals (Ag and Au) and the results were the 
same. In an effort to explain these results the 
violent-boiling theory of fragmentation was postulated (16). 
It was thought that perhaps the turbulence associated with 
the violent boiling of either the nucleate or transition 
regimes was responsible for initiating breakup. Then, 


based on Westwater's (21) demonstration that under 
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equilibrium conditions nucleate boiling occurs below the 
fluid critical temperature, fluids with a low critical 
temperature would tend to quench the molten material in the 
quiet film region allowing ole ee However, those 
with a high critical temperature may incur violent boiling 
prior to solidification resulting in fragmentation. This 
seems compatible with the observed results in water and 
sodium. It must be noted that this mode of fragmentation 
has not been observed in the case of uo. in sodium (17). 

In an isothermal investigation conducted by Ivins (22) 
consisting of mercury dropping into room temperature water, 
a Weber Number effect was observed. For values above a 
certain critical contact velocity fragmentation occurred. 
The fragmentation was considered to be a hydrodynamic 
phenomenon which occurred as a result of inertial forces 
Overcoming the surface tension of the drop. It was shown 
by Hinze (23), ina further study of the hydrodynamic mode 
of fragmentation that a critical value of the Weber Number 
results in fragments that are stable in nature. The 
attempt was to validate a method of fragmentation for coolants 
which exceed their critical temperature, but the results have 
not proven universally true. 

Cho (15) attempted to show the dependence of the extent 
of fragmentation on the initial temperature of the hot 
material and the coolant temperature for different molten 


materials dropped into water. The results demonstrated that 
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there is no consistent pattern as to the extent of 
fragmentation when altering the hot material temperature. 
It was concluded, however, that decreasing the coolant 
temperature results in an increase in fragmentation. 
Qualitative discussions of similar experiments (molten 
metal-water) reported by Brauer (24) and Flory (25), 
supporting the ideas of enhanced fragmentation with 
increasing Weber Number, and decreasing pool temperature, 
have yielded additional observations. No violent boiling 
occurred around the hot material and fragmentation took 
place almost immediately following contact. It was 
observed (25) that fragmentation resulted from an outward 
burst of the metal and that the residue was spongy in 
appearance. In the specific case of aluminum, many of the 
drops swelled into hollow thin-shelled "bubbles", some 
of which later burst. Further, chemically reducing the 
molten material's surface tension increased the tendency 
for fragmentation while similarly increasing the viscosity 
of the water reduced fragmentation. The resulting 
interpretation led the investigators to believe that 
encapsulation of the coolant in the hot material via 
Surface instabilities (Helmholtz instabilities) was 
responsible for the observed fragmentation mechanism. 
Alternative encapsulation mechanisms have also been 


proposed (19, 26, 27). 
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Amblard et al. (19) observed that in UO./water 
experiments no fragmentation took place. This was 
subsequently confirmed by Witte et al. (28). In an 
attempt to extend their results to the case of UO./sodium, 
Amblard et al. proposed that either an encapsulation 
mechanism or thermal stresses developed in the UO,, or a 
combination of both may be responsable for fragmentation in 
this system. 

In contrast to the above, some large-scale experiments 
consisting of dropping kilogram quantities of UO. into a 
pool of sodium were conducted by Johnson et al. (50). The 
UO. was heated to 5500°-5800°F via a thermite-type reaction. 
These experiments were similar to the small-scale experiments 
conducted at ANL (15) and the results were similar as well. 
Extensive fragmentation occurred but the interactions were 
mild in nature, and large pressure pulses of low energy 
content were observed. 

In a follow-on effort aimed at more precisely 
identifying the actual mechanism leading to fragmentation, 
stevens et al. (30) conducted experiments on transient film 
and transition boiling around spheres in water. It was 
observed that when the sphere temperature was lowered 
below a point where stable film boiling could be sustained, 
a very rapid collapse of the film layer took place. This 
rapid film collapse was suggested as a feasible method of 


fragmentation inducement. 
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Experiments consisting of dropping molten tin into 
water as reported by Board et al. (31) have revealed yet 
another set of observations. Various perturbations consisting 
of varying the tin temperature, varying the water temperature, 
changing the system pressure and introducing external 
mechanical impulses were carried out. The results led to 
the conclusion that fragmentation was triggered by the 
development of unstable film boiling between the two liquids. 
The main cause of dispersion was determined to result from 
vapor collapse which possibly generated a small coolant jet 
that penetrated the molten material, dispersed, and resulted 
in fragmentation. 

In an extension of his initial experiments, Board (31) 
investigated the effect of rapid vapor-blanket collapse. 

Once the molten material was in the coolant a rapid 
pressure increase was initiated by rupturing a diaphragm, 
thus collapsing the vapor blanket surrounding the molten 
material. An explosion was generated following the blanket 
collapse. These results were obtained using both tin and 
aluminum as the hot material and they seemed to confirm 

his earlier hypothesis. 

Other molten metal/water experiments conducted by 
Witte et al. (3la) using low melting point materials 
(mercury, lead, zinc, bismuth, tin and aluminum) appear to 
further substantiate Board's hypothesis. The results lead 


to a rejection of earlier theories, and resulted in the 
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conclusion that fragmentation iS a response to an external 
stimulus while the sample is still molten. It was further 
concluded that vapor film collapse was a highly likely, 
indirect initiating mechanism, while the direct mechanism 
for fragmentation was seen as being one of pressure 
differences, surface tension changes, or a thermal shock 
phenomenon. The results, however, were not sufficient to 


detail the exact cause of fragmentation. 


2.2.4 Coolant injection experiments 

Small-scale coolant injection experiments have been 
carried out at ANL utilizing various combinations of 
materials. Armstrong and Cho (15) injected small amounts 
of water at room temperature into a variety of molten 
materials. Both above and subsurface injections were 
Carried out, and in some cases explosive interactions were 
recorded. Anderson and Armstrong (32) noted that the 
existence of a gas layer between the two fluids seemed to 
prevent development of an explosion (similar to observations 
in Ref. 29 on the inverse experiment). 

Other experiments involving the injection of molten 
metal jets into water were reported in Ref. 29. The jet 
Surface temperature had a major impact on the extent of 
fragmentation. Additionally, when the jet was blanketed 
by a vapor film the initiation of fragmentation could be 


completely suppressed. Prior to the fragmentation explosion 





28 


a bulge in the jet was noticed (at the point of subsequent 
fragmentation), Similar to the phenomenon reported earlier 
by Flory (28). 

Small-scale coolant injection experiments have also 
been carried out utilizing sodium into uO.. The results 
are similar to those noted for water, as in some cases 
explosions developed, and a comparable time delay was 
observed. In this particular case it was suggested that 
Fragmentation of the coolant, followed by mixing with the 
molten material, resulted in the explosive interaction. 

An alternative explanation offered by Fauske (33) suggests 
that spontaneous nucleation of superheated sodium may have 
been the cause. He does not, however, exclude the 
possibility of an alternative mode of coolant fragmentation 


prior to Spontaneous nucleation. 


2.2.5 Shock tube experiments 

Another of the early efforts to ascertain the mechanism 
driving the high pressure levels encountered during a 
reactor accident led to a series of experiments designed 
SO as to impart a Significant thermal transient to the 
coolant via a heat source (34, 35). The results demonstrated 
that rapid heating of water from a flat surface and thereby 
producing vaporization yielded relatively low pressures. 
‘These pressures were even lower if non-condensable gases 


were present. Subsequently, the experimental methodology 
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was reviewed and a new type of experiment was devised: 
that of utilizing a shock tube. Shock tube experiments 
consist of supporting a column of liquid coolant above the 
hot material via a diaphragm, then puncturing or removing 
the diaphragm and allowing the liquid column to impact on 
the hot material. Wright (36) conducted shock tube 
experiments with water as the coolant and alternatively 
solid and molten aluminum, silver and UO. powder as the 
hot material. When solid materials were used, only low 
pressures were developed. However, when impacting molten 
materials, very high (nearly 6,000 psi) pressures were 
obtained, as well as dispersion (fragmentation) of the 
material. The resulting pressure was directly related to 
the molten material temperature, and the duration of 
pressure build-up was on the order of 0.5 msec. Hillary 
and Darby (1) later verified Wright's results and showed 
that non-condensable gas in the interaction zone (under 
low pressure) can stifle the explosive interaction 
altogether. 

Holtbecker et al. (1) utilized a sodium column ina 
shock tube impacting on alumina at 2670° K. Their results 
demonstrate repeated small-magnitude pressures too low to 


be destructive. 


2.2.6 In-pile testing 
Extensive fuel pin failure experiments involving molten 


fuel-coolant interactions have been performed as a part of 





30 


thermal reactor safety studies (37, 38, 39). Recently a 
parallel effort concerning in-pile experiments related to 
LMFBR systems safety has been taking place (40-47), and a 
summary of recent tests can be found in Ref. 48. The bulk 
of the LMFBR fuel pin failure studies have been carried out 
in the Transient Reactor Test Facility (TREAT). Data with 
respect to the threshold and mode of fuel-element failure, 
movement of fuel during transient heating, swelling, melting, 
and post failure effects of fuel-coolant interactions has 
been recorded. 

In an attempt to derive information relating to the 
pressures and the mechanical Sense generated in molten 
fuel-sodium interactions, ANL set up the stagnant sodium 
piston-autoclave tests (S-series). These tests provide the 
most Significant results relating to the fuel-coolant 
interaction in an LMFBR as a result of an over-power 
transient. Here the energy input yielding fuel melting and 
pin failure is actually supplied by a nuclear power pulse. 
The resultant fuel-coolant interaction involves fragmentation 
and pressure pulse generation. Yet, the energy conversion 
(nuclear to mechanical) is several orders of magnitude less 
than the maximum thermodynamic limit (49). Analytical 
interpretation of these experiments by Epstein and Cho (47) 
shows that the molten fuel was the only fluid capable of 
performing expansion work and that the liquid sodium acted 


aS an energy-dissipative fluid rather than a working fluid. 





Silk 


It appears likely that the initial pressure pulses and 
mechanical work measured results from fuel vapor pressure 
and bond-gas pressurization without involving any process 
of rapid sodium vaporization (48). Similar results have 
been obtained in the E and H series tests which use flowing 
sodium, and L and R series tests which simulate a loss of 
flow. That is, in all cases only mild interactions with 

a low energy conversion ratio have been observed. Further, 
the results of the UO.-sodium interactions experiments 
indicate that the energy conversion ratios are less than 
for the UO,-water case, for a comparable nuclear energy 
mow: (38, 39). 

The overall outcome of the TREAT experiments, although 
limited in scope, have provided valuable information and 
have advanced the understanding of fuel-failure mechanisms. 
It has been demonstrated, however, that large~scale vapor 
explosions are not a necessary outcome of a reactor 
transient. This remains to be verified analytically and/or 


experimentally. 


2.3 Theoretical studies 

In attempting theoretical investigations of the problem, 
two basic approaches have been undertaken. The first 
consists of attempting to parametrically model the overall 
thermal interaction from either a thermodynamic point of 


view, estimating the maximum work available, or froma 





32 


transient, rate-limited viewpoint. The second approach 
breaks the problem down to a subsystem level and involves 
studying only the possible mechanisms which lead to 
fragmentation of the fuel. In all cases it has been the 
intent to enable one to ascertain the energy conversion 
rates and subsequent mechanical work generated, the effect 
on the reactor system dynamics, and specifically the 


potential for large-scale vapor explosions. 


2.3.1 Thermodynamic models (parametric) 

The first attempt at determining the destructive work 
potential of a molten fuel-sodium interaction was carried 
out by Hickes and Menzies (49). It was a limiting case 
calculation which resulted in an upper bound on the 
expansion work available from the coolant. Initially upon 
contact the fuel and coolant were assumed to be in thermal 
equilibrium, sharing the thermal energy which had been 
previously stored in the core fuel. Then the resulting 
internal energy change which would be incurred by the fuel 
under an isentropic expansion to atmospheric pressure was 
taken to represent the work potential available. Later 
work has offered up modifications, perturbations and 
improvements to this original model (50, 51, 52, 53), with 
the most notable being that of Judd (52), which utilizes 
a more accurate equation of state for high-temperature 


sodium. 
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Although this type of thermodynamic approach provides 
an upper bound for the problem, it does not take into 
account the transient heat transfer and expulsion rates 
which actually occur, nor does it allow for the time 
dependence of the pressure rise or the mechanical work 
generation. As the problem is a dynamic one, the time 


dependence should be a key parameter in any realistic model. 


2.3.2 Rate-limited parametric models 

In an attempt to account for the transient nature of 
the fuel-sodium interaction, models which describe it in 
terms of time-dependent pressurization, voiding and 
rate-limited heat transfer have been developed. These 
models vary between those investigating a local interface 
phenomenon and those which assure that the fuel fragments 
and disperses within the coolant. 

Models of the limited interface type which have been 
proposed (54, 55) yield values on the opposite end of the 
spectrum from those of Hicks and Menzies (49). In these 
models, the heat transfer is assumed to stop following 
vaporization of sodium at the interface and the pressure 
work developed from a shock wave produced in the coolant 
is then evaluated. The speed at which vaporization takes 
place (less than 1 msec) accounts for the low work 
potentials developed. 

A model designed to incorporate the transient 


intertependence between the heat-transport and expansion 
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processes was developed by Padilla (56). He formulated 

the problem as consisting of spherical fuel particles 
uniformly distributed and in intimate contact with liquid 
sodium in an interaction region. It was further assumed 
that boiling was suppressed and that the liquid sodium was 
restrained by the remaining colder sodium from expanding 
when heated. The energy transfer to the sodium results in 
pressurization and subsequent acceleration of the surrounding 
liquid. Once the sodium pressure falls below the saturation 
pressure of the heated zone, the onset of boiling occurs. 

At this point the heat transfer to the sodium ceases and an 
isentropic expansion continues. During the heating process 
the only mechanism for pressure release results from the 
compressibility of the liquid plug above the interaction 
region. In this model, the effects of condensable and 
non-condensable gases are neglected, which could affect the 
peak pressures calculated. In addition, it may not be on 
the conservative side to assume heat transfer cutoff at the 
inception of boiling as it does not allow for energy 
transfer to the system during boiling. 

Building on Padilla's model, Cho et al. (57) at ANL 
developed a parametric model for the fuel-coolant 
interaction (ANL-FCI). A significant variation from the 
Original formulation was that of allowing for a rate-limited 
generation of fuel surface area, accounting for the finite 


fragmentation and mixing rate of the fuel in sodium. The 
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time-dependent area generation formula was not based on 

any physical picture of the process involved, and altering 
the generating time constant allows significant variation 
in the maximum pressures attainable (predicted). It was 
demonstrated that non-condensable gases reduced peak 
pressures as well as delaying its time of occurrence. To 
allow for the inclusion of non-condensable vapor blanketing 
of the fuel in thermal resistance, a later option has been 
added to the model (58). As a result of parametric studies 
using the ANL model (59), it has been demonstrated that 
noncoherence in the interaction leads to reduced mechanical 
work. However, this reduction is estimated to be small 

for realistic mixing times. 

Similar model development has taken place in other 
countries, and reviews of principal aspects have been 
presented by various authors (l, 2, 3, 60). The majority 
of these models have followed the above Cho-Padilla type 
formulation. The approach is to consider that particle 
dispersion has already taken place and that the fuel is 
contained ina finite volume of sodium prior to the 
interaction. Variations in the postulated heat transfer 
mechanisms, particle size and distribution, mixing modes, 
equations of state utilized, etc. contribute to the 
variations in predicted pressures. 

An initial attempt at the incorporation of a 


fragmentation model in the overall fuel-coolant interaction 
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calculation has been set forth by Grossgut et al. (61). 
Here, both solidification of the fuel and hydrodynamic 
impact forces are considered. Results appear to agree 
favorably with experiments. However, as pointed out by 
Caldarola (2), the number of parameters in the model may 
allow the generation of reasonable results just through a 
judicious choice of numerical values. 

The large number of optional parameters involved in 
any and all of the models generated thus far, coupled with 
the degree of uncertainty involved preclude blind acceptance 
of their results. They do, however, provide a convenient 
mechanism for sensitivity analysis of these various 
parameters. Two of the basic inputs which have a strong 
influence on the generated pressure and mechanical work 
are the contact mode and the fragmentation phenomenon, 
and until these can be determined with a reasonable degree 
of accuracy, the predictions of parametric models will be 


in doubt. 


2.3.3 Models for specific applications 

Models of this type are developed for a contact mode 
determined by the specific casualty sequence under 
consideration. A typical model of this type is that 
based on a jet-type ejection of molten fuel through an 
orifice in a failed fuel rod (62). Two others are the 


multichannel core dissembly code (FISFAX) (63), and ANL's 





oy 


SAS/FCI (64) model. The latter is primarily based on the 
ANL-FCI parametric model discussed previously. These models 
are strongly dependent on the underlying assumptions for the 
specific postulated accident, and their results are 


severely restricted in applicability. 


2.3.4 Fragmentation theory 

An excellent review of fragmentation theories/models 
to date has been compiled by Cronnenberg and Grolmes (65). 
As a result, only a generalized look into these models 
will be given here, with only significant details noted. 
As fragmentation and the subsequent mixing, contact, etc. 
of the fuel and coolant are major determining factors in 
evaluating the amount of mechanical work developed ina 
fuel-coolant interaction it is of great importance to 
understand these phenomena. As can be seen from the 
guantity and variety of theories proposed thus far it is 
evident that this understanding is largely lacking at 
present. 

It is convenient for organizational, comparative and 
discussion purposes to divide the existing models into 
three categories as per Ref. 65. These classes consist 
of: those models associated with hydrodynamic effects 
between the molten material and coolant independent of 
thermal conditions, the molten droplet models which assume 


that the liquid material is readily deformable and 
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undergoes fragmentation due to pressure forces initiated 
in the coolant, and the solid shell theories which consider 
thermal cooling effects and surface solidification. 

One type of hydrodynamic fragmentation model is that 
based on impact fragmentation as postulated by Hinze (23). 
This model, as previously cited, was based on a Weber 
Number dependence and the fact that a critical Weber 
Number corresponded to the onset of the fragmentation 
phenomenon. Another of the postulated hydrodynamic models 
is that based on Helmholtz surface instabilities (25, 66, 
67). Wave-like ripples which encapsulate coolant leading 
to breakup are postulated to form on the hot material's 
surface. These shear-induced waves depend on relative 
velocities of the materials and their magnitude and 
frequency depend on the material viscosities. A serious 
contradiction results as fragmentation has been shown to 
occur when only low or no velocities of consequence exist. 

The most widely postulated models are those based on 
a molten droplet form for the hot material. One mode of 
fragmentation falling into this group is that as suggested 
by Swift and Baker (68) involving vapor bubble growth and 
collapse. Considerable effort has been given to this mode 
of fragmentation, and variations thereon (7, 15, 69). In 
the dynamic bubble collapse mechanism set forth by Caldarola 
and Kastenberg (70) the destructive energy is considered 


as derived from the impingement of microjets generated via 
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bubble collapse. This formulation was a result of 
evaluation of photographic experiments (71). The jet 
impingement is descirbed as generating an elastic wave in 
the molten material which then generates acoustic energy 
resulting in fragmentation. 

The results of calculations based on these molten 
drop models, with sodium as the coolant, indicates that the 
maximum thermodynamic work votential is more than that 
necessary to account for the observed fragmentation. 
However, when estimating the actual amount of energy 
converted it is found insufficient to result in fragmentation. 
Thus, as fragmentation is known to occur, a more 
definitive estimate of the percent of energy conversion 
1S required. 

Another of the molten drop models consists of the 
generation of an acoustic pressure oulse in the coolant 
(69). The pressure pulse is a result of thermal expansion 
of sodium at the fuel-coolant interface as a result of 
heat transfer. However, again the maximum work 
potential is smaller than that required for fragmentation. 

Still another possibility is that of the previously 
noted concept of spontaneous nucleation of the coolant 
following liquid-liquid contact as considered by Fauske 
(72). Whereas this mechanism is plausible for some 
metal-water interactions, the UO.-sodium contact 


2 


temperature is lower than the spontaneous nucleation 
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temperature of sodium. It has been pointed out, however, 
that for large amounts of uo. Surrounding a small sodium 
droplet the requirements for spontaneous nucleation can be 
met (73). 

A final molten drop mechanism is that based on internal 
acoustic cavitation in the hot material (14). In this 
case surface boiling generates internal pressure waves 
which then produce internal cavitation leading to breakup. 

Epstein (75) has proposed a violent gas release 
mechanism from a metastable supersaturated solution as a 
fragmentation mechanism. This model is in a category of 
1ts own, and is highly dependent on the solubility of 
gases in the molten material. As of now the mechanical 
work potential of this model has not been calculated. A 
Significant fact is that such a mechanism is capable of 
accounting for the sponge-like appearance of some 
fragmented materials. 

In contrast to the molten drop models, recent efforts 
to take into account the possibility of solidification of 
the hot material have resulted in solid shell theories of 
fragmentation (76-79). Initially, Hsiao et al. considered 
the pressurization of a solidifying sphere of aluminum in 
an infinite water cooling medium. It was concluded that 
the tangential stresses will exceed the radial stresses 
and that the maximum surface tangential stresses occur as 


solidification begins, which should allow for immediate 
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rupture. Ref. 77 extended Hsiao's results to the case of 
vO. and sodium with similar results. The cornerstone of 
these models is the solidification kinetics of the hot 
material. For the case of UO,-sodium (78, 79), the time 
for molecular reordering was found to be short compared 

to the time constant for heat transfer. Fragmentation is 
assumed to result from coolant penetration of fissures 
developed in the outer frozen shell or from boiling within 
these cracks. The fact that the surface stresses developed 
in the solidifying shell are found to exceed the yield 
strength of UO. by several orders of magnitude led to the 
conclusion that fracture should occur. The actual 
fracture mechanisms, as well as the energy conversion 


rates, and the particle size and distribution rates have 


not as yet been determined. 





CHAPTER THREE 
THERMAL STRESS GENERATION 


foe introduction 

The solid shell theory of fragmentation, under 
consideration here, appears to be consistent with the 
limited experimental observations (i.e., Armstrong, 
Ref. 20) for the case of U0.-Na. The work conducted to 
date (75-76), although limited in scope, has shown it to 
be both feasible and relevant. As pointed up in the 
literature review, the fracture mechanisms, time scales, 
particle size distribution, etc. have not, as yet, been 
determined. In order to make use of this fragmentation 
theory in an overall Fuel-Coolant Interaction model, these 
quantities must be known. Thus, the ability to determine 
them in a straightforward analytic manner is highly 
desireable. Prior to making an attempt at extending the 
theory, a brief review of the solid-shell fragmentation 
theory is necessary to obtain a first-order estimate of the 
desired quantites. This existing theory will then provide 
the foundation for fracture calculations and time scale 


estimations. 
42 
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The analysis consisted of making use of Adams's (82) 
approximate solution to the heat conduction equation for a 
constant surface temperature to determine the liquid-solid 
interface velocity with time. The constant temperature 
boundary condition was determined from the contact interface 
temperature formulation of Ref. 83. Utilizing molten UO. 
initially at its melting point (2800°C) and 200°C sodium, 
an asymptotic solid-liguid interface velocity of 2 cm/sec 
was found to exist for times greater than 1 msec (Figure 
3.2). This result was based strictly on heat transfer 
considerations. 

Secondly, classical crystallization theory was used 
to determine the rate of crystal growth in the melt (see, 
for example, Ref. 78). Based on the maximum solid-liquid 
interfacial surface tension the homogenous nucleation 
temperature was determined to be 700°C below the melting 
point of UO.. A more realistic estimate of the solid-liquid 
interfacial energy was made and this value was then 
revised to be a few hundred degrees. If any nucleation 
Sites are activated prior to this temperature, heterogenous 
nucleation would be initiated and solidification would 
“commence immediately. Nonetheless, for the UO.-Na 
temperatures under consideration, the contact interface 
temperature is well below even the lower limit homogenous 
nucleation temperature and solidification is sure to occur. 


An investigation of the various temperatures associated 
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figure 3.1--Model for solidification induced stresses. 
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with all realistic LMFBR accident conditions revealed that 
this would also be the case. The net rate of crystallization 
for the 1218°K contact interface (the same as that used in 
the heat transfer calculations) was found to be 71 cm/sec 
(Figure 3.3). When this value was compared to the mean 
solidification velocity of approximately 2 cm/sec found 
from heat transfer conditions, it was concluded that the 
uo. solidification is limited by the heat transfer process 
and not the time for molecular reordering. Consequently, 
it can be assumed that the solidification begins immediately 
upon contact of the molten uO. with sodium. 

Armed with the fact that solidification of UO, will 
occur, and the knowledge that Hsiao (76) had demonstrated 
the feasibility of thermal stress fragmentation, an 
evaluation of the pressurization and thermal stresses in 
the UO.-Na system was carried out by Cronenberg et al. (77). 
In this case, the effectsof temperature dependence of the 
mechanical properties, compressibility of the molten liquid 
core, and variation of the surface heat transfer conditions 
were investigated. In conducting parametric evaluations 
of the validity of various assumptions, both Al-water and 
UO.-Na systems were analyzed for comparative purposes. 

The solution for the temperature-dependent thermoelastic 
stress problem of Ref. 84 was used for rotationally symmetric 
stress and strain. It was further noted that for the case 


of constant properties and a Poisson's ratio less than 0.5, 
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@2  sOlid shell fragmentation theory 
The initial attempt at analytically modeling the 


fragmentation of solidifying hot material in a coolant was 
that of Hsiao (76). His investigation was centered around 
determining the pressurization of a solidifying sphere of 
molten material in an infinite cooling medium (in his case, 
aluminum and water, respectively). In order to attack the 
heat-conduction problem (a partial differential boundary 
value problem with a nonlinear surface boundary condition) 
the steady-state approximation of London and Seban (80) was 
used. To gain the desired approximate solution, the 
following assumptions were imposed (76): 
1. The molten material was taken at the uniform 
melting temperature, initially. 
2. Thermal properties were taken as constant. 
3. The unit surface conductance, h, was taken as 
constant. 
4. The temperature of the liquid phase was 
considered uniform and constant. 
5. The hot material was considered to have a 
discrete melting temperature. 
6. The density of both phases was taken as the 
same. 
In addition, the molten material was taken as a small 
Sphere, suddenly immersed in the coolant; solidification 


was assumed to begin instantaneously, and the compressibility 
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of the liquid core was neglected. The geometry of the 
Situation is as illustrated in Figure 3.1. 

Based on the foregoing assumptions, the steady-state 
approximation was then employed to obtain an evaluation of 


the rate of solidification front movement and the 


temperature profile in the solidified shell. Then, 


calculating the induced pressure and the thermal gradient, 
the pressure and thermal stresses were evaluated. The 
maximum total stress was found to occur ina tangential 
direction at the surface of the drop, almost immediately 
after solidification began. It was concluded that if the 
temprature was such as to yield stress levels which exceed 
the ultimate tensile stress of the shell, failure would 
Secur. 

The stress calculations were based on elasticity 
theory, and in light of the observations of Ref. 81 
concerning the high strain rate and rapid temperature 
change, modeling the material behavior as elastic was 
considered adequate. 

Before extending this type of model to the UO,-Na 
system, Cronenberg et al. (78, 79) investigated the 
validity of the assumption of solidification initiation 
upon contact. The analysis was based on a calculation of 
the crystallization kinetics in an effort to determine if 
the time for molecular reordering to form solid from the 


melt is short compared to the heat transfer time constant. 
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the pertinent equations reduced to the familiar form (85) 
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where T(r,t) is the temperature difference between the 
present and initially unstressed states. 

In order to determine the thermal stresses, both the 
temperature and the velocity of the solidification front 
need to be evaluated. For constant thermophysical properties, 
the Fourier heat conduction equation in spherical geometry 


was formulated as 


2 
x mee (r,t) im 


or 


elit, t) | 2 or (r,t) (3.4) 
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with a' the thermal diffusivity, and the initial condition 
being 


Au, © ) = T (3755) 
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The ligquid-solid interface conditions are given as 


Jap ete) or, Galt) (3.6) 


k ox r=a dt 


for an energy balance, and 
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T(a,t) = Th (3.7) 
Bor continuity. 
The external surface boundary condition is 


Oy Ui) 


=k or R 


= h(T. (t) = T,) (328) 


where the coefficient of heat transfer, h, is dependent on 
the mode of heat transfer. 

In contrast to the convective boundary condition, 
the previously mentioned constant contact interface 
temperature condition could exist if perfect wetting occurs 
(a highly likely situation). In this case, the interface 
temperature is given as 

i Usa Cia) ie To (k./vegjert 2 
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where H and C denote the hot and cold materials respectively. 
Incorporation of the erf A’ term results from the movement of 
the solidification front in the vo. and its value is 


determined from the transcendental equation (77) 








C T 2 Kv de 
Sa = ie + erf } (3.10) 
LY vO. Kava, 





for vo, liquid initially at its melting point. The values 


for bulk UO. at its melting point in 200°C sodium are given 


2 
as: A = 0.94, erf »X = 0.82, and the temperature of the 
interface is 945°C (77). The contact interface temperature 


Will remain approximately constant for as long as the uO, and 


Na remain in contact and the entire uo. sphere has not solidified. 





a2 


The number of approximate solutions is extremely 
limited, especially for spherical geometry. In the study 
under discussion here, two separate approaches were taken. 
First, a numerical solution of the equations was utilized 
in an effort to include transient effects. Second, an 
analytic solution based on the constant surface temperature 
condition was employed. The constant surface temperature 
solution utilized the approximate solution of Adams (82) 


with the temperature in the solidified shell given as 


Z 2 
4 L R-r| 2 da date wae pe be R=5 
at) = ca? a iE | “ det a ada a Fg Bz] 





(3.11) 
Rorlat yay 
a9 «6 ‘ye 3Rr 
with the freezing rate given as: 
? -_ 
20, 5 Ws Ty) 
Ga _ (oil) 
at 
Me sa2aylis{i+42 2-3) 
3 aL m iE 








The position of the solid-liquid interface at any time is 
found by integrating the above equation. 

Of major interest in the initial heat transfer 
analysis was the determination of the proper cooling 
surface boundary condition and its range of applicability. 
A comparison was made between the constant surface 
temperature condition and a convective nucleate boiling 
heat transfer coefficient. The selection of the nucleate 
boiling coefficient for use was based on the prediction of 
Henry's correlation (86) which required a temperature of 


9500°C for Na film boiling. Figure 3.4 illustrates the 


JY 
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effect of changing the surface boundary condition on the 
liguid-solid interface position. As stated (77), the 
solidification of vO. in Na is limited primarily by its 
low conductivity and not by the heat transport process, 
and the choice of the cooling surface boundary condition 
is not extremely critical for h = We lieste Seat ame 
especially over short times. 

Several other parametric evaluations were carried 
out by Cronenberg, as well (77). A comparison of the 
transient numerical solution with the steady-state 
approximation of London and Seban (80) used by Hsiao (76) 
was made (Figure 3.5). The results indicated that 
neglecting thermal inertia, as the steady-state approximation 
does, results in deviation of the solution. It was further 
found that using temperature independent properties (with 
mean values of a and E) or varying Poisson's ratio from 0.3 
to its upper limit of 0.5 had little effect. Similarly, 
the limited compressibility of the molten liquid core only 
produced a small decrease in the internal pressurization. 

As was found in Ref. 76, the total tangential stress 
at the surface was the dominant factor in the U0.-Na case, 
as well. Figure 3.6 illustrates the thermal, pressurization, 
and total stress components at the outer surface as a 
function of time. The total surface tangential stress is 
Shown to go through a maximum at approximately 50 msec 


after quenching, a value in line with the observed breakup 
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time of uo. in Na observed experimentally (20). This 

point coupled with the fact that the total stress at the 
surface exceeds the yield stress of UO. (5000 psi at 

2000 C) by a Significant amount led to the conclusion 

that the thermal stress theory of fragmentation in the case 


of UO, in sodium is indeed plausible and it deserves further 


2 


consideration. 





CHAPTER FOUR 


FRACTURE THEORY 


4.1 Introduction 

The intent here is to briefly review some basic 
theories of fracture of ideal solids, the fundamentals of 
the fracture mechanics approach to fracture analysis, and 
the alternative surface energy criterion for fracture in 
ceramics. As all work on the fuel-coolant interaction to 
date has been based primarily on heat transfer/thermodynamic 
type analyses, a review of these basics was considered 
essential prior to attacking the modeling concept considered 
here. Although this treatment is somewhat Superficial, it 
is believed that the concepts and terminology will provide 
sufficient understanding where application of the 


principles of fracture is required. 


4.2 Crack propagation in homogenous materials 


In an ideal, homogenous, isotropic, elastic solid, 
the tensile strength is on the order of its molecular 
cohesion force. While the strengths of small whiskers, or 
carefully produced fibers, are found to approach this value, 
the experimental strengths of bulk material are generally 


only 1/1000 - 1/100 of the ideal value. The reason for 
58 
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this discrepancy, as explained by Griffith (87), is the 
presence of small flaws or cracks which result in severe 
stress concentrations in even the unstressed state. fThus, 
even though the average stress in the solid may be 
relatively low, the theoretical fracture strength can be 
locally exceeded. 

When attacking the problem from an energy standpoint, 
it can be seen that when a material is stretched, the 
interatomic Spacing of planes of atoms is increased. This 
increase in distance iS proportional to the stored elastic 
energy. It is this potential energy, as ascertained by 
Griffith (87), that provides the surface energy required to 
form new surfaces during fracture. For most materials, 
other than the most brittle, it has since been found (88) 
that a Significant amount of the apparent effective surface 
energy 1S absorbed in plastic flow and the formation of 
branch cracks. Furthermore, when fracture propagates at 
high speed, energy is lost in the production of sonic 
waves (88). 

Griffith's analysis was for a plane flat plate of two 
dimensions with a stationary central crack, elliptical in 
Shape. He allowed the crack to become very sharp at its 
ends. Then, assuming Hooke's Law valid, he extended 
Inglis' analysis (87) to calculate the critical nominal 
stress. Inglis described the locally magnified stress at 


the crack tip as: 
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= 2 
ee 20 /E (4.1) 
where SP xe stress at crack tip 

0 = applied nominal stress 

cs = half crack length 

0 = radius of curvature at the crack tip 


The radius of a sharp crack tip is difficult to measure, 

and so ey cannot easily be evaluated by this formula. 
Griffith therefore extended the analysis based on energy 
considerations. From the stress distribution around the 
elliptical crack, he determined the stored elastic strain 
energy. Then, equating this stored energy to the energy 
required by the newly formed surfaces, he found the critical 


stress to cause crack propagation (for plane strain 


conditions): 
where 
i Critical applied nominal stress cauSing 
fracture 
E = Young's modulus of elasticity 


c = crack half length 
y = fracture surface energy 
v = Poisson's ratio 


Under conditions of plane stress, this equation becomes: 


ee ey (4.3) 
S 1C 
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In the above, plane stress is defined as a state of 
stress which is described in every volume element by a set 
of principle stresses, one of which is zero. Generally, a 
plane stress condition exists in thin sections in which 
the stress normal to the surface is nearly zero. On the 
other hand, plane strain is defined as a state of stress 
which results in a zero strain along a specific direction. 

Prior to the initiation of crack growth, in the 
Griffith theory, the system strain energy increases 
directly with the applied stress. Once the growth begins, 
however, the change in strain energy becomes a function of 
the magnitude of both the crack opening kinetic energy 
and the energy required to form new crack surfaces (surface 
energy). The kinetic energy here results from the fact 
that as the ends of the crack move forward, the material 
at the sides of the crack moves apart with a finite 
velocity. Thus, a kinetic energy can be associated with 
this movement of material near the end of the crack. This 
kinetic energy, initially zero, increases as the crack 
grows. This relation can be graphically illustrated as 
shown in Figure 4.1, as a function of the half crack 
length. 

The locus of the stress and strain at the onset of 
crack instability can be illustrated as shown in Figure 
4.2. Then, where the stress-strain curve of the material 


intersects this locus, fracture should occur. Typically, 





62 


KINETIC ENERGY 


Co 
CRACK LENGTH 


Figure 4.1--Kinetic energy versus half crack length (89). 
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Figure 4.2--Stress-strain relation based on 
Griffith criterion 
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the point of intersection occurs in the region where the 
slope of the locus is positive. The possible discrepancy 
between observed fracture stress values and those predicted 
by the Griffith equation (4.2) is related to the angle 
between the stress-strain curve of the sample and the locus 
at the point of intersection. With increasing crack 
length, the value of this angle approaches ninety degrees, 
where the Griffith criterion is satisfied. The 

implication is that, due to the difference in the strain-energy 
stored prior to fracture, a small crack would have a higher 
starting velocity than a large one. 

Further, the kinetic energy, fracture surface energy 
and the strain energy can all be determined graphically 
from the fracture locus (89). Consider a sample containing 
a crack of length 2c which is extended to fracture, with 
the ultimate stress maintained during crack propagation. 

If the path OAB shown in Figure 4.3 is followed, the breakdown 
of the work into its various energy components, at any time 
following the onset of crack propagation, can be determined 
mig. 4.3). 

Sack (90) and Sneddon (91) extended the Griffith 
criterion to a two-dimensional crack (disk shaped) ina 
three-dimensional solid. The basic assumption underlying 
this analysis is that the tensile strength of the brittle 
material in one direction is unaffected by the stresses 


at right angles to it. Then the highest allowable stress 
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Figure 4.3--Graphical determination of 


fracture surface energy and 
kinetic energy (89). 
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| 
| 
| normal to the plane of the crack is given by: 


om — tt .,.. (4.4) 
2c(l-v~) 
where 
nominal stress applied perpendicular to 


the plane of the crack to cause fracture 


E = Young's Modulus of elasticity 
y = fracture surface energy 

v = Poisson's ratio 

c = radius of the disk (crack) 


Depending on the value of Poisson's ratio for the material, 
. from equation (4.4) differs from Griffith's solution by a 


maetor Of 1.57 to 1.8 (90, 91). 


4.3 Crack propagation velocity and branching 


The crack propagation velocity is dependent upon 
various factors, including: the material, the imposed stress 
pattern, and the amount of stored strain energy at the time 
of fracture initiation. There is a velocity limit, referred 
to as the terminal velocity, and it iS a material property. 
When the available energy exceeds that necessary to propogate 
the crack at its terminal velocity, the excess will appear 
in other forms or it will cause crack oscillation or 
branching. 

In comparing theoretical and experimental values (92) 
of the terminal crack velocity in brittle materials, it 


has been found that the governing factor is the kinetic 
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energy supply. The extent to which the kinetic energy is 
supplied to various parts of the body determines the terminal 


velocity. The terminal crack velocity is then given as 


(92): 
v = 0.38 /—+. /1 - = (4.5) 
9) (Cha) 
V = terminal crack velocity 
E = modulus of elasticity 
o = density of the material 
a = Original length of the crack 
c = extended length of the crack 
v = Poisson's ratio 


It has been found that the terminal crack velocity 
is independent of the surface energy, and that the stress 
state in the material surrounding the crack is approximately 
equal to that in the static case (92, 93, 94). The 
foregoing solution yields a terminal velocity of approximately 
38% of the velocity of sound in the material: /E/p. 

The relationship between crack velocity and the stress 
required to sustain it can be seen in Figure 4.4: the 
stress decreases with increasing speed (88). Therefore, 
once a crack is initiated, it may continue to propagate 
even though the stress level decreases. It is possible 
then to have catastrophic propagation even under a rapidly 
decreasing load (i.e., decreasing from a maximum). Since 


the crack tends to be self-sustaining at high velocities, 
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Figure 4.4--Tensile stress to sustain crack velocity 
versus propagation velocity (92). 
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it is reasonable to assume that branching must occur 
prior to the occurrence of this condition. 

It 1s possible to obtain a solution for a two-dimensional 
crack in a three-dimensional medium by evaluating the 
stresses and displacements in a semi-infinite solid (95). 
In this case, it has been found that the terminal crack 
velocity is again independent of the surface energy (96). 
The conclusion was that a two-dimensional crack moves in 
one plane, with no branching, at a constant velocity. 

When the Griffith criterion is satisfied by the 
stress condition in a sample, the sample is in an unstable 
state and crack propagation begins with an infinitesimal 
increase in the stress. Practically speaking, the stress 
would then continuously increase until the sample breaks. 
The actual fracture stress value depends on the loading 
rate and the Specimen geometry. Figure 4.5 illustrates 
the relationship between the crack length and the elapsed 
time following crack initiation. The time required to 
obtain a given crack length decreases as the applied 
stress increases. The reason for this behavior results 
from the different initial velocities of crack propagation 
for different stress levels. However, the terminal velocity 
is the same for all levels; this is illustrated in Figure 
4.6. Figure 4.6 also shows that for lower stress levels, 
a longer crack is required prior to attaining the terminal 


velocity. 
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Figure 4.6--Crack propagation velocity versus crack 
length (89). 
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It has been demonstrated from experiments on glass 
(97, 98) that the terminal velocity of all cracks is the 
same, unless the crack slows down or stops as a result of 
branching. The crack propagation velocity, however, is not 
uniform across the specimen width (99). The velocity of 
each individual fracture is influenced by such factors as 
the degree of load relaxation during the dynamic phase of 
the process. 

In a specimen subject to bending, the velocity 
history 1s somewhat different. The crack initiates and 
propagates through a zone of transition, but as the 
terminal velocity is approached, the specimen's neutral 
axis cannot shift fast enough as the crack approaches. The 
resulting velocity versus depth is shown in Figure 4.7: 
the velocity decreases with depth under flexure. The 
decrease results from the fact that the crack has passed 


beyond the original neutral axis location. 


4.4 Basic fracture mechanics 

As most recent trends in crack propagation analysis 
have centered around the local stress and strain relations 
at the crack tip via what is known as fracture mechanics, 
a brief review of this concept is considered essential. 
This type of analysis is based on the belief that crack 
growth will occur when certain critical conditions are 


reached at the crack tip (100). Essentially, the 
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Figure 4.7--Fracture velocity versus depth for a glass 
bar in bending (100). 
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application of fracture mechanics depends on knowing the 
material behavior under various crack loading conditions 
from experimental results. Then, comparing this behavior 
to the applied conditions at the crack tip of any crack 
a determination as to whether the critical conditions are 
reached can be made. 
4.4.1 Brittle fracture (linear 
elastic fracture mechanics) 

As previously related, Griffith's original work (87) 
for fracture of an ideal brittle material was based on 
an energy balance criterion. This work was later extended 
by Irwin (101) and the Griffith equations were rewritten 
in terms of a crack extension force, g, (required for 


crack movement) as: 


g = = (plane stress) (4.6) 


neo” 2 
aa (Oe 


(plane strain) (4.7) 
where 
c = the half crack length 
o = the applied stress 
v = Poisson's ratio 
This crack extension force must be greater than a critical 
crack extension force for crack propagation. 
Griffith introduced cracks of various lengths and by 


recording the loads at failure he was able to determine 


that the product of the original crack length and the load 
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at failure was essentially a constant for a given 

material. Irwin (101) then related the crack extension 
force to a stress intensity factor, K, which took the 

place of Griffith's experimental constant in the stress 
equations describing the stress field around the crack tip. 
Physically, K may be regarded as the intensity of the load 
transmittal through the crack tip region resulting from 

the introduction of the crack, and it therefore represents 
the magnitude of the stress field in this region. Irwin's 


relations were given as (101): 


Z 
ge — (plane stress) (4.8) 
K? 2 
—— (l - v*) (plane strain) (4.9) 


The symbol K_ is normally used to represent gv/1c 
(the critical value) and it is called the fracture 
toughness of the material. For a given material, all 
combinations of crack geometry and loading at incipient 
fracture yeild the same value of Ky making it a constant 
for that material. The fracture toughness, or critical 
stress intensity factor as it is sometimes called, Kye can 
be a function of temperature or loading rate (99), and it 
represents the critical value of the stress field in the 
immediate vicinity of a sharp crack tip as illustrated in 
Figure 4.8. 

The general form of the equation representing the 


Crack tip stress field in terms of K, when the load or the 





76 


"(€OT) PTaetF ssoays dtQ yoezD. au ut squauodwods ssaizj4s ayy pue 
yoezrd ayy Jo abhpa Hhutpest sy worFJ paeanseouwl saqzeutTprACOD--g°y, aanhtg 


WIVYD SHL JO 3903 ONIGVSI 





77 


crack size is below the threshold of instability is (101): 





K 

oe = — Be (8) (4.10) 
where 

mk = X,Y,Z 

n = failure mode (I,II,IITI) 

EA = a function of the angle 

K = constant, independent of r 

fe) = stress in directions i,k 

ie = distance from crack tip 

v) = angle 


The fracture modes are denoted as (I) wedge opening, 
(II) forward shear, and (III) parallel shear as depicted 
in Figure 4.9. 

The factor K is analytically determined, and it 
varies as a function of the crack configuration and the 
method of load application. However, once the critical 
value of K for a material is determined, for a given 
combination of crack configuration and loading, it remains 
essentially unchanged for all other combinations of crack 
and loading conditions. Hence, asS previously pointed out, 
the critical value of K is referred to as a material 
constant. An excellent compilation of the analysis of 
cracks and the determination of K for numerous cases and 


geometries can be found in Ref. 103. 
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When analyzing stresses and fracture properties under 
real loading conditions, it is necessary to take into 
account all types of loadings. For different loadings 
resulting in the same failure mode (stress intensities 
caused by pressure, temperature, etc.) the total stress 
intensity factor for the system is the sum of the single 
stress intensities of the same mode (102). As an example, 
for mode I this implies: 


KL = K, (1) ot K, (2) oc 


This procedure and the above formulas are strictly 
applicable only in the linear elastic field. 
4.4.2 Elastic-plastic fracture 

(plasticity effects) 

In an effort to determine the influence of yielding 
at the crack tip, other models and methods of attack were 
developed. The procedures representative of such an 
analysis are given by Wells (104). Here, the plastic zone 
around the crack tip is described by: 


2 





oS K 5 (4.11) 
Y 2110 
y 
where 
a = plastic zone size 
o., = yield stress 


Then from the displacement equations, the crack opening 


displacement is found to be: 
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ea noe 
§ - Eo Eo (4.12) 
y 
where 
6 = crack opening displacement 


In this case, a critical crack opening displacement 
represents fracture. 


In a more general form, the plastic zone is given 





as (103): 
Ne 
ry = z (4.13) 
y 
where 
— (pla tress) 
| ah plane stre 
OL — 
i 
67 (plane strain) 


The plasticity at the crack tip causes some 
redistribution of the stresses to maintain equilibrium. 
Thus, the full width of the plastic zone is taken as twice 
the above result. These results are somewhat in error, as 
they neglect work hardening, large strains, or other 
influences. In some cases 7 is added to the crack length 
to obtain an “effective crack size" followed by conventional 
linear elastic analysis. This method is not extremely 
accurate, but has proven useful when trends are desired. 

However, it has been pointed out (103) that small 
amounts of plasticity or other non-linearities at the 


crack tip do not seriously further disturb the load 





} 
31 
redistribution around the crack tip, hence the relevance 
of the K values remains. Therefore, where small amounts 
b Be non-linearity are deeply embedded within the stress 
field, they do not have a significant impact, and a 


combination of K K and K represents a unique stress 


io £1 JGIGAE 
field at the crack tip for smali-scale yielding. MThus, if 





| 
yielding is limited to a small enough scale, the 


linear-elastic theory is both reasonable and adequate. 
When plasticity effects are present to the extent 
such that general yielding is in evidence, the elastic-plastic 
_ analysis may not be completely valid. Several basic 
approaches in the area of non-plane strain behavior are 
being developed. One such criterion being studied at 
Westinghouse (104), originally developed by Rice (106), 


is based on an energy line integral, J, defined as: 


7 2 |. wdy - T 2 ds (4.14) 
where 

T = traction vector 

u = displacement vector 


w = strain energy density 
The application of J is similar to that of K in the 
linear-elastic analysis. As yet, this method has not been 
developed to the point where it can be generally applied, 
and presently elastic-plastic analysis methods are used 
aS approximations in many cases where extensive plasticity 


1s found to exist. 
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4.5 Fracture of ceramics 

Ceramics, compared to other materials, are extremely 
brittle. They are thus prone to catastrophic failure, 
particularly under conditions of shock. As a result, 
ceramics have come under separate study. One such study 
is that of Ref. 107, which provides a unified approach in 
determining the fracture of ceramics. 

As noted before, the application of the Griffith 
criterion yields values of the fracture stress too small 
when plasticity has an effect. Since for most materials, 
the high shear stresses near the crack tip exceed the 
plastic flow stress, allowing some stress relief, the 
applied stress must be increased for propagation to continue. 
To account for these plasticity effects, the alternative 
energy criterion has been presented, as opposed to the K 
(stress intensity) type of approach of fracture mechanics. 
According to this::approach (107-109), the effects of 
plastic flow can be incorporated ina surface energy for 
fracture initiation, Yue so that 


5 yee 


On = eo (plane stress conditions) (4.15) 


In theoretical evaluations, allowances are made for all 
energy losses due to various forms of work (plastic flow, 
cleavage steps, crack branching, etc.) but crack blunting 
is not accounted for. Here For ok it is usually not a 


sufficient condition for crack propagation. It would be 
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sufficient, however, if an empirical Ys which accounted 


£ (except 


at high temperatures) and the effects of crack blunting 


for crack blunting is used. In ceramics, os eG 


are likely to be small so that the measured Y,; May 
approach the theoretical value (110). 

If the size of inherent flaws can be determined, the 
above equation can be used to calculate the stress to 
extend these flaws. If the values obtained are close to 
the observed fracture stress, it is likely that the 
fracture is a direct result of extending inherent flaws 
(brittle fracture). If this is not the case, fracture is 
most likely initiated by plastic flow. Here, the fracture 
stress is related to the stress to nucleate and/or 
propagate fresh flaws. Then the fracture control 
mechanism requires an analysis of the plastic flow stress. 
The former case, however, normally holds true for ceramic 
materials. 

The alternative surface energy criterion is related 
to the K (stress intensity) concept in that at the :point 
of catastrophic crack propagation the following holds: 

Ke = ar = 2Ey;- Utilization of 0 in the case of ceramics 
is preferred in some cases because it can be directly 

related to the energy absorbing process occurring during 
crack propagation (110). In metals where blunting by plastic 


flow takes place, K is more useful. 
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Moretical relationships do not exist between Ve and 
grain size, porosity, temperature, or purity, and there is 
not enough data to obtain empirical relationships. 
Therefore, the utilization of this method is, at present, 
somewhat limited. Variation of VR with crack length has 
been evaluated for a number of materials (110, 111). At 
low to moderate temperatures Vea is found to be independent 
of crack length for edge cracks greater than a grain 
diameter in length. 

In estimating the theoretical fracture stress for 
brittle failure via extension of internal flaws, the relevant 
flaw size must be known. This flaw size, in ceramics, is 
related to the microstructure of the material through 
pores, voids, grains, etc. The largest flaw is controlling, 
if it is atomically sharp. When pores are larger than 
grains, they can act as flaws, themselves. In most cases, 
the inherent flaw size can be related to either pore or 
grain sizes. 

Normally, the stress to extend a small flaw through a 
grain is much lower than the fracture stress. Therefore, 
in these cases, the inherent flaw existing in a ceramic 
can be approximated by the pore size plus the grain size 
(111). Then with the value of Ue known, the fracture 
stress can be evaluated, or vice versa, through Equation 


Be 14. 





CHAPTER FIVE 
FRACTURE ANALYSIS OF THE UO.-SODIUM SYSTEM 


Sel introduction 

In order to determine whether or not the solidifying 
shell of UO, will fracture, various methods of making such 
a calculation were investigated. In all cases, the 
analysis was necessarily quasi-static. After evaluating 
fracture stress, energy balance, and fracture mechanics 
type methods, an approach based on fracture mechanics 
was determined to be the most appropriate. This method, 
based on the knowledge of the stress intensity factor, Kis 


and the fracture toughness, K will then allow 


idee 
determination of whether the shell will fracture froma 
given flaw, under the influence of the prevailing stress 
field. A description of the foundation of this type of 
fracture analysis was given in the previous chapter. In 
this case, due to the high strain rates and the brittle 


nature of the UO the linear-elastic approach was 


2* 
considered adequate for first-order calculations. 

In order to obtain immediate and inexpensive 
evaluations of the possibility of brittle fracture, the 


application of analytical methods is most desirable. 


These analytic methods can be checked by means of finite 
85 





86 
element computer techniques, if the application warrents. 
The result is that there are three basic ways of attacking 
the problem; the first two being analytic and the third 
the finite element method. First, a one-shot approach 


based on existing analytic models assumes that once K, 


exceeds Koc catastrophic crack propagation takes place. 


The second, also analytic, is based on varying the crack 


length across the shell, ensuring that Ky exceeds Koc all 


the way across (a progresSive method). The third method 
is the same as the second, but it uses finite element 


techniques to evaluate K. rather than existing analytic 


it 


models. The first method was selected for the initial 


investigation. 


5.2 Analytical evaluation of Ky 


As shown previously, it is possible to compute the 
stresses caused by the internal pressure and by the 
temperature distribution. The components of these stresses 
perpendicular to the plane of the crack are needed for the 


computation of the stress intensity factors K (resulting 


It 


from the thermal gradient) and K (resulting from the 


Ip 
induced pressure). Brittle fracture will then occur when 


the sum of these stress intensity factors exceeds the 


fracture toughness, K 


IC 
The thermal stress distribution across the shell can, 


for a given time step, be linearly approximated as shown 


in Figure 5.1. The linearized thermal stress distribution 





37 


TANGENTIAL STRESS 





———— approximate 
actual 





Figure 5.1--Linear avvroximation of thermal 
and pressure stress components. 
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Figure 5.2--Linear approximation of total 
stress. 
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can then be combined with the pressure stress distribution 
(approximately constant) at the same time step. The 
combined stress distribution can then be completely 
described by stress values at two radial positions since 


the distribution is assumed linear (Figure 5.2). 


Cmax °D 4 “thermal max (DoH) 
Omax Z Omin 
Cmean 2 seg) 


In order to calculate the stress intensity factors 
in terms of the applied stresses, an analytical expression 
is then required. The linear thermal stress distribution 
is similar to that resulting from an applied moment. Thus, 
the two stress fields can be represented as the superposition 
of a uniform stress upon that caused by a bending moment. 
An analytic expression for just such a case has been 
developed (111) to allow for the solution of an arbitrary 


linear stress field. The total K. is defined as (with a 


18 
representing crack depth in this case) 
a a 
2S = = mee 
Rit | nant (3 a omean' ere) 


where g(a/w) = k(a/w) - £f(a/w). 

This expression, along with the appropriate 
correction factors, is summarized in Figure 5.3. 

An alternate method based on the analytical models 
given in Ref. (103) can also be utilized. In this case 


the values of Kip and Kit are calculated separately and 





CORRECTIONS 





O of O2 03 04 O8 O6 
ae 
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K; = (tra)"'? | Onay f (a/w) +Gmean 9(a/w) | 





Figure 5.3--Stress intensity factors for a single 
edge-notched strip in combined stretching 
and bending (102). 
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added together to yield the total K Again, the two 


r° 
stress fields (pressure and thermal) are modeled as a 
uniform stress field and a linear stress field equivalent 


to that resulting from a comparable moment, respectively. 


These models are illustrated in Figures 5.4 and 5.5. 


®.3 Comparative evaluation of the modeling technigues 


In an effort to compare the two alternative analytical 
models, a sample calculation was carried out to evaluate 
the predicted Ky values. The stress field values in the 
solidified shell used in the comparison were those for a 
time of .044 sec as reported in Ref. (78) (see Figure 5.6), 
as this time is close to that at which the maximum 
tangential stress level is reported to occur. The 
linearized approximation used in the analysis is also 
shown in Figure 5.6. Further, an approximate flaw size 
er 75 x pon: cm was used, representing a grain size of 
25 um and a pore size of 50 um as obtained from reference 
(111). This value is an estimation based on the fact that 
in polycrystalline UO. the stress to propagate a pore to 
fracture is accurately predicted by taking a flaw size 
equal to the pore size plus the grain size (lll). The 
comparative calculations for each model can be seen in 
Figure 5.7. As shown, the difference in the results is 
within the accuracy of the calculations. Therefore, the 
first method was arbitrarily selected for use in succeeding 


evaluations of Ky values. 
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LINEARIZED STRESS APPROXIMATION 
(COMPRESSIBLE LIQUID) 


Figure 5.6--Actual and linearized approximation 
for the stress distribution in the 
shell. 
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Method 1 
= 4.7 x 10° psi = 1.55 x 10° psi 
max ° P Cmean : : Ee 
a/w = .167 > g(a/w) = .40 and Cay wee lL 
” 1/2 
K, = (ma) [onax f (a/w) + 006,79 (a/w) I 
K, = ao X oe psivin 
Method 2 
3 5 S 5 
op = imo x LO psi Ce 3.2 x 10° psi (max) 
Pressure 
Kip = Opyma F (a/b) 
(a/b) = .167 + (1 - a/b)?! “PF (a/b) = .49 
F(a/b) = 1.509 
Kip = 2.18 x 104 psivin 
Thermal 
Ki. = ° F (a/b) (a/b) = .167 
B72 = 
(1 - a/b) F(a/b) = .8 > F(a/b) = 1.062 
a, = 3-273 x 10° psivin 
° a = 4 — 
y Ky = Kip + Ki = 5.45 x 10 psivin 


Figure 5.7--Analytic evaluation of stress intensity factor, 


Kir via two alternate methods. 
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5.4 Determination of Koc 


Through the analysis of extensive amounts of data on 
the deformation and fracture of oxide fuels (113, 114) the 
brittle fracture stress of UO. has come to be represented 
as 


mee 6,029 + 3.3337 (5.4) 


where T is in °C and o, iS in pSi. The input data behind 


ie 


equation (5.4) over the range of interest was established 


Via three-point bending tests of 24 x 6 x 5 mm uO. 


Specimens (111). Note that UO, specimens were prepared 
by powder pressing and sintering methods and may not be 


characteristic of UO. shells obtained from solidification 


2 
of molten, irradiated UO. However, for this material it 
is expected that the grain size and porosity characteristics 
are Similar, and that the resulting fracture characteristics 
Will not vary significantly from the ones used here. It 
is then possible, utilizing a nominal flaw size of 75 um, 


to determine the appropriate values of K with respect to 


ne 


temperature. The simplest "best" estimate of Kc 1s 


obtained through the use of the analytical model illustrated 
in Figure 5.8. Application of this model, in conjunction 
With the fracture stress found from equation (5.4), yields 


the plot of K versus temperature as shown in Figure 5.9. 


ine 
To engage in first order fracture calculations, a value 


corresponding to a mid-range temperature is selected. In 
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this case a value of Koc = 2.407 x no psolyane vac 


approximately 2000°C) was chosen. 


5.5 Evaluation of the stresses in the solidified shell 

In an effort to obtain the stress states and their 
appropriate values for a variety of times, the fixed surface 
temperature condition, using mean property values (77) was 
employed for the conditions reported in Chapter 3, i.e., 


UO, at its melting point, Na at 200°C and an interface 


2 
temperature of 945°C. These temperature conditions were 
maintained for the entire UO./Na analysis in this chapter. 
Based on Adams (82) solution for the temperature profile 
and solidification front velocity (equations 11 and 12 of 
Chapter 3), equations 1-3 of Chapter 3 were used to 
evaluate the induced pressure and stresses. The induced 
pressure, in this case, is based on the attempted thermal 
contraction of the solid shell with the limited 
compressibility of the molten uO. taken into account. The 
density change upon solidification is neglected. The 
resulting values of the solidification front position with 
time, induced pressure with time and the external 
tangential stress with time can be seen in Figures 5.10- 
me 12. 

The values obtained via the above method differ 


somewhat from those reported in Ref. 77. In the present 


case, the overall stress levels are found to be higher 
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Figure 5.10--Liquid-solid interface position with fixed 
surface temperature. 
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(although of the same order of magnitude) and the peaking 
phenomena experienced by the external tangential stress 

in Ref. 77 is nothing more than a minor perturbation here. 
This latter result came from the fact that the maximum 
tangential pressure stress Occurs almost instaneously as, 
initially, the pressure rises faster than the shell 
thickness increases. The values of the thermal stresses 
developed in the shell are higher than those of Ref. 77 
mor all times, including the initial value at t = Ore This 
initial value is directly related to the values of a and 
E used in the stress equations and the values used here 
are consistent with those reported in Ref. 77. Thus, it 
is not clear wherein lies the difference. Hence, the 
values of the induced stresses obtained in the present 
Study are taken as those to be used in the fracture 


analysis, tO maintain consistency. 


5.6 Variation of the surface heat transfer coefficient 
In an effort to determine the effect of various 
amounts of entrapped fission gas surrounding the UO. 
particles in the sodium coolant on the generated stresses, 
an analysis based on varying the surface heat transfer 
coefficient was also carried out. Taking constant, mean 
values of the thermophysical properties, the governing 
equation and conditions are given by equations 4-8 of 


Chapter 3. In conducting the analysis, the steady-state 





E03 


solution of London and Seban (80) was used for the heat 
transfer and solidification calculations in a manner similar 
to that used in Ref. 76. This method was considered 
adequate as only a comparative evaluation over relatively 
short times was desired. The same assumptions as enumerated 
earlier with respect to this method are applicable, and 
the geometry of the model is still that of a molten vo, 
sphere suddenly immersed in an infinite sodium pool (Figure 
el). 

The steady~state approximation, strictly valid when 
the latent heat of fusion is significantly greater than the 
internal energy change in the shell, maintains the shape 
of the temperature profile as solidification proceeds. 
The governing equations for this case, replacing 4-8 of 


Chapter 3, become: 








2 
d_(rt(r)] _ 4 forea <1 <.-K (5.5) 
2 
dr 
where T(r) = T'(r) - TO subject to the boundary conditions 
~k wat = hT(R) at r=R (5.6) 
= r=R 
T(a) = T =T' - T at r = a(t) C5257) 
m m C 
alt x ) d (a) = 
k = _ oL ag = at xr a(t) (5.8) 


mma a(Q) = R. 


The solution of the above is 
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aha Ps [2] [2 ap 4 [ar 7 1} 
"BS 


stax ~ 2] 


in nondimensional form. Equation (17) yields the temperature 








profile within the solid shell. Substitution of equation 


(17) into equation (16) results as 


in oe =a eee (Se0) 
[afl 
[Fler - 1) + 2] 
where 
r= a/R 
and 
— 2 
t = kT t/LR 
m 


Then, integrating subject to the initial condition 
fe=) (a/R) = 1.0 at t = 0 
yields the equation for the solidification front position 


eee | (e> = 1) 


Bi re 


et ae 2 eo 
3 + — =o = t (5201) 


A plot of r vs. t for various values of h demonstrates 
the relative sensitivity of the solidification rate to the 
coefficient of heat transfer, and thus, the amount of 
Surrounding fission gases (Figure 5.13). The constant 
Surface temperature results are included as well for 
comparative purposes. It can be seen from Figure 5.13 
that as the heat transfer coefficient is increased, the 
rate of solidification increases. Further, it is noted 


that, as concluded by Cronenberg et al. (77), the fixed 
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surface temperature condition results in values close to 
(although slightly faster than) those corresponding to a 


~ 


nucleate boiling coefficient (41.34 Cal/sec ae 8 ze 
mo,000 Btu/hr ft-°F), for the same coolant temperature 


with UO. at its melting point. 


2 
The second part of the comparison, and of equal 
importance, consists of determining the effect of varying 
the heat transfer coefficient on the stress levels at the 
external surface (point of maximum stress). Again, only 
a comparative evaluation is required, so the steady-state 

temperature profile was considered adequate. 

Substituting equation (5.9) for the temperature into 
equation (3.3) and integrating yields the pressure 
developed in the molten core. The result, written in 
non-dimensional form, is 


{(3)-(8) J-2(l)-la) le 


= i a een SD) 
E(T -T ) 3 3 . 
m c a ie a Za a 
[1+ 2 [gr- 2] |{a+2 ty-ay {3 uy z [1-(3 } 


and the tangential pressure stress component of equation 











2 
R 


m.2) is 


J ae 3 
eerowier + KR) ieee 


O — 
Pt de (ae = R?) 


Similarly, substituting equation (3.7) into the thermal 
stress component of equation (3.2 yields the tangential 


thermal stress as 














(5.14) 


alE-tlé-21 
i falar] 


Utilizing the values of the solidification front 


al 
2 


position with respect to time over the range of interest 
(as shown in Figure 5.14, which is an expanded view of 
Figure 5.13), the total surface tangential stress was 
calculated, and plotted as seen in Figure 5.15. Not until 
the value of h falls below approximately 10 cal/sec cm? °C 
does any significant variation occur. The interesting fact 
is that even though the pressure developed in the molten 
core decreases with decreasing heat transfer coefficients, 
the solidification rate is retarded, creating a thinner 
Shell which yields an increased pressure stress. The thermal 
Stress, on the other hand, decreases with decreasing values 
Of h. However, as h decreases, the pressure stress begins 
to dominate the thermal stress and rises rapidly enough to 
Cause the total stress to increase over the initial value 


based on a nucleate boiling h. As was seen in Figure 5.13, 
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the fixed surface temperature condition results in faster 
solidification. As a result, this boundary condition 
should provide the highest total stress level anticipated. 
The values of the heat transfer coefficient were 
originally selected to reflect a maximum (41.34 Cal/sec 
em°C) roughly corresponding to that for nucleate boiling 
of sodium and a minimum (.1 Cal/sec cm*°c) below that which 
would correspond to film boiling. The manner in which these 
compare to approximate values for trapped fission gases 
is illustrated in Figure 5.17. Figure 5.17 is based on 
the approximation tnat 


k 


~ s = (5,125) 
where 

k = thermal conductivity 

§6 = film thickness 


The values for thermal conductivity employed were selected 
with the value for helium as an upper bound and a value 
1710 of that for helium for a lower bound. The value of 
k for helium was taken at a mid-range film temperature 
from Figure 5.16 which was based on (115) 
k = 0.002418(991 + 0.678((T/1.8 - 1200)) (50) 
for 
Hog eR < TF < 6GOO0O°R 
The overall range of film thicknesses which roughly 


Corresponds to the values of h used is from 4 x tome to 





3.0 


-3 


2.0 


K = 0.002418 (991 + 0.678 (( T/1.8 — 1200)) 


1600 °R < T < 6000°R 
Ref. 86 


K, THERMAL CONDUCTIVITY ( Cal/sec-cm-°C) x 10 





2 3 


TEMPERATURE (10° °C) 


Figure 5.16--Thermal conductivity of helium. 
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> xX om cm. This range is considered broad enough to 


adequately represent the values anticipated in the actual 


UO.-Na system. 


2 
From this study, it can be concluded that the 


entrapment of fission gases does not appear to reduce the 
stress levels in the solidifying shell and in fact can 


cause them to actually increase. 


5./ One-shot fracture 
The one-shot fracture method assumes that once Koo 
determined through the use of an analytic model, exceeds 


the critical value, K catastrophic fracture ensues. 


ne. 


Thus, taking the stresses described in the last sections 
and the analytic model illustrated in Figure 5.3, the value 


of Ky at a given time can be computed and compared to Kic 


to determine if fracture occurs. 
Taking as a representative example the stress field 


in the shell for the UO.-Na system with the interface 


temperature fixed, at t = 5.76 x ome sec (see Figure 5.18), 


a value of K, = 7.67 x 107 psi Yin is found to exist for 


an assumed length of 75 um. Since this value of KS is 


greater than K crack propagation, according to this 


es 
theory, will occur. Further, since the magnitude of the 


Stresses does not change significantly over the first 100+ 


millisec, K, will always exceed K over this range for 


a IC 


the same size crack or larger. Therefore, any flaw 
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Figure 5.18--Total tangential stress distribution across 


the shell. 
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generated in the shell corresponding to a size of 75 um or 


more would prove catastrophic. 


5.8 Progressive crack growth evaluation 


In an effort to verify the one-shot approach and to 
ascertain whether or not crack arrest may occur, an 


analytical evaluation of K. with the crack progressing 


If 
across the shell was undertaken. The model employed in 
this case was based on that of Ref. 116 for calculating 

the thermal stress intensity factor based on the 
Superposition method. This model was initially derived 

for calculating the stress intensity factor for a plate 
with one face subject to sudden cooling, and with a flaw 
at or near the surface being cooled. The thermal stresses, 


J, (x) | in the crack free plate (determined from the 


at 
momentary temperature distribution) are used in conjunction 
with the known solution for a crack in an infinite plate 
loaded at its surfaces by a pair of wedge forces to obtain 
Ky via integration over the crack. The analytic solution 


for this case is given as 


Cc 
Ky pele LZ 2ve | 0, (x) | oe (Sess) 
C YT Wee = x" 


where c is the crack depth. The y direction is parallel 
to the plate surface, while the x direction is parallel 
mo the penetrating crack. The factor 1.12 describes the 


influence of the stress free surface where the crack 





116 


starts. The effect of the finite plate width is not 
considered in this formula. 

In this case, as with the one-shot approach, the 
possibility of brittle fracture is dependent on the actual 


crack tip stress intensity factor, K exceeding the 


ied 
critical stress intensity factor of the material, Ric: 
As Kio may be locally temperature dependent, and as KT 


is a function of geometry (crack deovth) and the time 
dependent thermally induced stresses, the basic criterion 
for rapid crack extension can be expressed as (116) 

K, (90 


meee] K(t,c) > Kk. (0) (spel) 


Using a guasi-static analysis, the variation of Ky with 
crack depth can be determined for each time step. 
Extending this model to the case at hand, that of a 
sphere rather than a plate, should yield a reasonable 
first-order approximation. Thus, assuming the crack to 
be long with respect to its depth, the region of the shell 
containing the crack can be approximated as a rectangular 
flat plate. The stress intensity factor can then be 


determined from the stress field in the undisturbed shell 


via the superposition method, and equation (4) becomes 





Cc o, (r) 
K, (c) = alee) 2 2vc See ee se aR - xr) (5.19) 
t VT ; az - (R- x)? 


As a result of the connectivity of the sphere, a 


Sertain amount of constraint prohibiting the complete 
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redistribution of the stresses around the crack is 
inherent. Therefore, an additional moment must be 
incorporated in an effort to allow for the conservation 
of the moment across the cracked section. The result will 
be an added effective stress distribution acting on the 
crack, causing the value of K, to increase with crack 
depth. Some methods have been proposed to incorporate 
this effect (117) but the validity of such methods is 
somewhat suspect where, as in this case, the tensile 
stress 1s a strong function of position. As a result, 
such a correction factor has been omitted in this work. 
The tangential stress field in the solidified shell, 
for various positions of the solidification front (l.e., 
time steps), can be seen in Figure 5.19. The Ky values 
corresponding to these particular times can then be 
determined and illustrated as shown in Figure 5.20 (note 
that the crack position in the shell is non-dimensionalized 
for plotting convenience). The significant factor 
illustrated by Figure 5.20 is that the stress intensity 
factor exceeds the critical value, Ric = 2.407 x 107 
psivin, almost instantly and it remains so thereafter. 
The time to fracture is then related to the time response 
of the thermal effects, and as soon as solidification 
begins and stresses are generated, fracture will ensue. 


This seems to verify the initial outcome of the one-shot 


approach discussed earlier. Note that the reason for the 
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rapid drop in the value of K_ as the crack is made deeper 


if 
for a given time step is a resuit of the strong radial 
dependence of the stress field. 

If, upon penetration of the shell by the initially 
activated flaws, the internal pressure is relieved, it 
may be possible to have the value of KS for remaining flaws 
fall to zero within the shell. If this were the case the 
arrest of subsequently generated flaws is possible. fThus, 
the next step was to investigate this last hypothesis. 
Figure 5.21 was established using only the value of the 
thermal stresses through the shell, i.e., assuming that 
the pressure had been relieved. As crack propagation only 
results if the crack tip stress intensity factor is greater 
than the critical value, the decrease in the Ky curves for 
the case of thermal loading only would seem to imply that 
cracks of growing depth would be arrested as KS falls 
be low Kio: This is only true if the crack exhibits stable 
growth, where no excess kinetic energy is available to 
drive the crack further. The exact point of arrest, 
directly influenced by dynamic effects, does not lend itself 
to simple analysis. The actual stability of the crack 
extension is dependent on the initiation mechanism and 
energy conversion rates involved. 

@rack stability is difficult to accurately predict. 


Various studies have been carried out (118, 119, 120), 


and thermal stress crack stability has received some limited 
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attention (119). As yet, however, no theories are 
available for crack stability and propagation in 
transient non-uniform thermal stress fields as exist in 
the present model. In general, if upon propagation the 
elastic energy release rate exceeds the surface energy 
required to form new crack surface, the remainder is 
converted into kinetic energy. The crack, once 
propagated, will continue to grow until all the kinetic 
energy and further released strain energy has been 
converted into surface energy. 

The rapid rise of the stress level in the shell 
for the case at hand, along with the magnitude of those 
stresses leads to an assumption of unstable growth. A 
qualitative assessment of Refs. 118-129 tends to support 
this theory. Thus, even if pressure relief occurs with 
propagation of the initial flaws, subsequently initiated 


flaws are anticipated to be catastrophic as well. 


5.9 Effect of grain size and porosity 
Saebrittle fracture strength 


Various factors can have an influence on the fracture 
Strength, and in turn the Kic value, of a material. In 
the case of brittle ceramics, such as U0., the two most 
important parameters in this category are grain size and 
porosity. The effects of these factors are not known in 


detail, but general experimental trends have resulted in 


@ealitative assessments of their influence. 
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One attempt at establishing an empirical relationship 
for the combined effect of grain size and porosity is that 
of Knudsen et al. (121), which is given as 


o = Kerams (5.20) 


where S is the fracture strength, G represents the mean 
grain size (um), P represents the fraction of porosity, 
and K, a and b are all empirical constants. This 
formulation is based on overall experimental trends, and 
it becomes coupled with a third influencing factor, 
temperature. Generally, an increase in temperature tends 
to strengthen specimens with a low porosity and a fine 
grain size, and to weaken those with a high porosity and 
a large grain size. Also, an increased sensitivity to 
changes in porosity and grain size has been noted at high 
temperature (122). 

In many cases, the effect of grain size on the 
maximum stress and/or the onset of plastic strain appears 
to be relatively small (123). This apparant contradiction 
Can be attributed to the fact that Knudsen et al. (121) 
were concerned with critical flaws on the order of the size 
of the grains, whereas the experimental observations of 
Refs. (111) and (123) were based on much larger critical 
flaws; thereby the grain size effect became less important. 
Furthermore, in the brittle fracture region, under 
consideration here (comparable to the low temperature 


region of testing) it is extension of the pre-existing, 
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inherent flaws which determines the fracture stress. Thus, 
it is the flaw size and not the grain size that becomes of 
major importance. 

The specific effects of porosity depend on such factors 
as: whether deformation is elastic and fracture brittle, 
fracture is preceeded by plastic flow, or extensive 
plasticity occurs. Further, there will be a critical 
strength dependence on pore morphology, and the operating 
conditions of temperature, stress, and/or strain rate. 

In the region of interest, that of brittle fracture, 
an increase in porosity yields a decrease in strength. 

The magnitude of this effect is determined by pore size, 
shape and distribution (124). Large pores allow for 

easier fracture initiation and crack propagation, and the 
stress to extend these flaws determines the fracture stress. 
This porosity effect on the fracture stress is also 

related through the influence on the modulus of elasticity. 

Another interesting effect results from the fact that 
at high strain rates, increased porosity yields brittle 
behavior (124). This is even true at moderately high 
temperatures. If the strain rate is decreased, however, 
plasticity may develop. As well, a density relationship 
with porosity has been shown to have an effect on the crack 
behavior of UO. (124). In the case of high density, a 
Single through crack develops, but with lower density 
Material, many part-through cracks were observed. A point 


Beenote is that not only could stress relief stop the 
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Gracks, but the onset of plasticity could allow crack 
blunting and subsequent arrest. 

As to the effect of grain size and porosity on the 
present analysis, the impact is not considered to be very 
Critical. As the strain rate is high and the fracture 
mode brittle, grain size effects should be minimal. 
Further, from the data of Ref. 111, the slight increase 
in fracture stress for a smaller grain size would offset 
the slightly smaller critical flaw size in the calculation 


®t K yielding essentially the same value. 


EC ’ 
A variation in the porosity could have a more 
Significant impact. The greatest effect would be in 
altering of the average pore size, which, being the 
Critical flaw size, alters the fracture stress. However, 
the direct dependence of Ric on both flaw size and fracture 
stress tends to dampen the variations as there is an 
inverse relationship in the brittle region (121). This 
fact, along with the knowledge that the induced K.. values 
Ic (requiring large variations in Koc @ 


make any noticeable difference) have led to the qualitative 


greatly exceed K 


conclusion that the Kic value employed herein is adequate. 


Meno Density change effect 


When vo, solidifies, there is a change in the 


Material density, and this perturbation on the problem is 


unaccounted for in previous sections of this report. In an 
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effort to determine the relative magnitude of the 
resulting effect and the consequences in light of the 
developing stresses, a parametric evaluation was conducted. 

The first estimate was based on determining the 
magnitude of the developing void which would form if the 
outside radius of the drop remains fixed. A further 
Simplification was made by assuming that the void developed 
uniformly in the liquid and had no effect on the 
progression of the solidification front. This may prove to 
be a poor approximation, as the void would be concentrated 
at the upper portion of the spherical drop (gravity/buoyancy 
effect) and thus inhibit solidification in this region. 
However, as the precise shape of the developing void is 
unknown, this approximation was considered adequate for 
initial estimates. 

For this first case of Ro (outside radius) fixed, 
the initial mass of the drop is given as 

4 3) 


= —— ae R 
m V TI a 0 


2 PV, 5 (S28) 


QR 
After a subsequent time step, the mass of developed solid 
is 


4 3 3 
— — <= Rae 
m 5 7 (Ro a 0. (5 ) 


and the volume of solid 


= 


V — _S 65.23) 
Ss ae 


Defining the void volume as 





eZ 


boid = Me 7 Ve a = (5.24) 


ama the void fraction as 


; “os dg Ws ~ Sh ~ Me 
void fraction = v = ae air ae (3225) 
fe) Oo 





substitution and manipulation leads to the expression 


{ 


Ps 


ee (5.26) 
Pe 


void fraction = /l - r 
( O 





Equating the original mass to the mass of the solid yields 
an estimate of the final internal radius for complete 


solidification with Ro fixed (assuming the void is maintained 





at the center). This expression is given as 
iG 8 
Pe 
a=Ril- >a (5.27) 
O 0 
S 
and for the current value of Ro = .36 cm, a = .2124 cm 


which yields a final void of approximately 20% by volume. 

The interesting aspect of this void formation 
phenomenon is that in its voresence, no internal 
pressurization will develop, and thus, no pressure stresses 
will be present in the shell. 

In an effort to bring the initial estimate closer to 
reality, and to gain more insight into the actual effects 
On the pressurization, the thermal contraction of the shell 
was brought into play. The change in the internal shell 


radius, due to thermal effects alone, is given as 


Aten 7 


ta] 


[oyy - Uo + O.,)] + aAT(r)a (By 233) 
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It is noted that the limited compressibility of the liquid 
would inhibit complete thermal shrinkage until a void 
actually began to form. The first estimate was to calculate 
the volume reduction due to thermal contraction, and then 
reducing the generated void by this amount, to determine 
the actual void production rate. The plot of void 
generation rate for both the case of fixed Ry and the 

case of considering thermal contraction can be seen in 
Figure 5.22. Note that two curves are given for the latter 
case; one is void fraction based on the original outside 
Roe and the other (perhaps the most pertinent) is based 

on the thermally contracted Ro: 

The plot of Figure 5.22 demonstrates the fact that 
there would be an initial delay time prior to any void 
formation. This results from the fact that the original 
shell contraction rate exceeds the rate at which the 
density change would cause a void to form. As a result, 
Maethe early stages of solidification, pressurization 
will occur. Once the void begins to form, the pressure 
would rapidly decrease and only thermal stresses would be 
present in the shell. The net result is that the initial 
shell stresses will be high enough to successfully 
propagate an existing flaw. Once the pressure is relieved, 


the situation would become the same as the non-pressurized 
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case investigated in the fracture section. The 

only difference is one of higher order, in that the 
Shell structure in the region of the void would no 
longer be symmetric, and there would be a temperature 
field perturbation. However, the fracture 
calculations conducted for the non-pressurized case 
Can be considered adequate for a conservative 


estimation. 


5.11 Constraint developed when Na 


becomes trapped in UO. 


As described in the review of proposed models, 
when the contact interface temperature is greater 
than the homogenous vapor nucleation temperature, 
spontaneous nucleation of the coolant is possible. 
This fact, in some cases, could yield a vapor 
explosion. Normally, as illustrated in Figure 5.23, 
for uo. in sodium, the contact temperature is well 
below the spontaneous nucleation temperature. However, 
as pointed out by Fauske (74), this mechanism 
might apply to the case where a small amount of sodium 
becomes trapped in molten UO, and is subsequently heated 
to its homogenous nucleation temperature. It is of 
interest, therefore, to investigate the inverse problem 


of Na trapped in vo. to determine the amount of constraint 
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supplied by the solidifying uo. shell when the sodium 
reaches the point of spontaneous nucleation. Thus, 

a crude, first-order estimate was carried out for just 
this purpose. 

The physical problem is postulated as illustrated 
in Figure 5.24. The internal sodium droplet is assumed 
spherical in shape and since it has a high thermal 
conductivity, its temperature is considered uniform, 
at an average value, for solidification calculations. 
As well, the external solidification front 
temperature is assumed to be constant at the vo. 
melting point, and the steady-state approximation of 
London and Seban is again utilized to evaluate the 
rate of advance of the solidification front. The 
method employed 1S a quasi-static analysis based on 
an energy balance. At each time step, the energy 
balance allows estimation of a new sodium temperature 
which can then be used to determine the solidification 
front movement to the next time step. Once the Na 
has exceeded its spontaneous nucleation temperature, 
the process is terminated. 

The energy change in the sodium can be approximated 


as 


AE... = Me AT - IV ,¢, (T, ~ Ti) Iya (5.29) 
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MOLTEN UO. 


Solidification Front Movement 





SOLID UO, 


LIQUID Na 


Figure 5.24--Liquid Na trapped in molten UO 
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Eguating the energy changes, substituting and manipulating 


leads to the following: 


3 _ - eae 
R ap. u = UT, _ Psolid@p 3 LB) Neher fo1d 


) ) 
2 1’ Na 4 Sind 
(5.32) 


The solidification equation, similar to equation (19) for 


the inverse problem now becomes (79) 





—- 1 il =n eee 
where 
Yr = a/R 
and 
t= -T kt/LR* 
a Na - 


Initially, to estimate a reasonable size and time step 
to use, a preliminary calculation was carried out. In 
this case a sodium drop is assumed to be instantly 
Submerged in a pool of molten UO,, and the sodium is 
assumed to have a negligible internal temperature 
gradient. Also, the effects of UO, Solidi rileaeloneane= 


neglected. For this case, the energy balance is written as 





ee 


An Bees) na - Ty oh) Ne dt (S24) 


which yields 




















Voc T t 
ml A inl - 
Na f Z 
OER © 
Na p eae TE 
t= = —N2 gn|-=—— (5.35) 
é 2 
where 
Te = surrounding uo. temperature 
a initial Na temperature 
T. = final Na temperature 


Taking a value of the sodium spontaneous nucleation 
temperature of approximately 2000°C (Figure 10) plots of 
time to spontaneous nucleation versus the surrounding fluid 
temperature for various Na drop radii and versus drop radii 
for various values of h can be constructed (see Figures 
9.25 and 5.26). Based on this initial investigation, 
values of R = .3 cm and h = 10 Cal/sec aC OE were 

selected for use in the first-order investigation 

including solidification. Then, employing equations (26) 
and (27), and a quasi-static analysis, Figure 2.27 can 

be developed. This figure shows the average Na temperature 
and the solidification front variation with time. Noting 
the time at which the spontaneous nucleation temperature 
(determination of this value will be subsequently 


explained) of sodium is exceeded (t = 1.25 msec), the 
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Figure 5.25--Time to spontaneous nucleation for Na versus 
temperature of surrounding fluid for varying 


drop size. 
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Figure 5.26--Time to spontaneous nucleation for Na versus 
drop size for varying surface heat transfer 


coefficients. 
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Mesition Of the solid-ligquid interface can be determined 
(a = .315 cm), and the uo. shell thickness is then 
defined as .015 cm. The value of the spontaneous 
nucleation temperature shown in Figure 5.27 is based 
on Figure 5.23 with the saturation temperature and 
induced pressure being determined as described in the 
following. 

The pressure induced in the sodium during the 
solidification of the uO. is computed neglecting the 


compressibility effects and the thermal expansion of the 


Na (offsetting effects). Similar to the development in 


Ref. 76 the stress effects on the radius of the solidifying 


shell can be given as 


aa ds _ 
Seeman be cet re * See)! F OpttyR 
Ar epee — vic. + 6..)1 
press Ee tr np tp 
taking ey Y oress yields 
60,,,E * 5 
P = a a ai a | ays jee > slits 
RX + 2a” - v(4a~ - R) R 


which, as before, reduces to (in nondimensional form) 
2 8) 
renee ellie st 1 
P 7 R| R R| |(R Ba 
aE (T —-T ) a Ve: 
mm 6|OCC a a le 
{(8) +2-»[4 a} ]} fs - 1} 


Utilizing mean constant values of the thermophysical 


vad 
R 








properties, equation (5.38) leads to the development of 


Figure 5.28, which snows the pressure history during 


(5.36) 


(55370 


Qe) 


(5.368) 
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solidification. The peaking phenomenon is a result of the 
heating of the trapped Na, gradually reducing the 
temperature difference across the shell as solidification 
proceeds. The decreasing temperature difference reduces 
the amount of thermal shrinkage and thereby decreases the 


induced pressure. The vapor pressure versus T. curve for 


at 
sodium is based on (125) 
wee Oy = _ 
log Eee = 7 Oia NG Give a9 2.5 5 (on. 59) 


Where T is in °K, and it can be seen in Figure 5.29. 
Pigures 5.28 and 5.29 can be used in conjunction with 
Figure 5.23 to determine the spontaneous nucleation 
temperature of sodium with time which is the line plotted 
in Figure 5.27 and previously described as defining the 
solidification limit. 

Finally, the stresses developed in the solidifying 
Shell can be evaluated to estimate whether or not the 
Shell would remain intact when the sodium reaches the 
spontaneous nucleation temperature. Thereby, a first-order 
estimate of the constraint at this time can be obtained. 
In order to estimate the temperature distribution in the 
shell, a quasi-static evaluation was carried out based on 
an isothermal surface boundary condition and constant 


conductivity in the shell with the temperature given as (84) 


af R a 
i -. R fr, a[2 - 2 - T.R{I - 2] (5.40) 
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Fiyure 5.27--Na temperature and solidification front 
position versus time for Na trapped in UO.. 
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PRESSURE vs TIME FOR 
Na TRAPPED IN UO, 


9 R = 03cm 
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Figure 5.28--Pressure induced in Na droplet 
trapped in solidifying vO... 
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Figure 5.29--Sodium vapor pressure as a function of Toat’ 
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Based on the above temperature distribution and the 
previously evaluated pressurization (Figure 5.28), the 


total tangential stress is then calculated from 





ee 3 2 
“- — — 
-_ = eo ee T(r,t)xr° lig Sax =< T(r,t)r- dr 
a il Vv 3 oe 3 3 
mia = KR )r ZY 
R R 
- = T(r,t)| + eRe a) 
2r {a ~ R ) 


which is equation (3.2) rewritten for the inverse problem. 
The results of the stress calculations are shown in 
Figure 5.30 (plotted for the maximum value which is at the 
internal surface). Note that previously, even when only 
considering thermal stresses (Figure 5.21) in the original 
problem, the Koc value of 2.407 x ne psi/in was 
substantially exceeded, i.e., by an order of magnitude. 
Even though the stresses found in this evaluation of the 
inverse problem are less than those yielding Figure 5.21 
(approximately 5 x on ps1), they are of the same order of 
Magnitude, and thus fracture is again anticipated as 
occurring. The time to spontaneous nucleation was 
estimated as 1.25 msec (Figure 5.21) for the case at hand. 
The fact that the stress is high enough to generate Ky 
values sufficient to cause fracture for all times less 
than this leads to the conclusion that fracture will be 
instantaneous (i.e., prior to spontaneous nucleation). 


Thus, no solid shell constraint would exist at the time of 


Spontaneous nucleation. 
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Pefore leaving this subject, it is of interest to 
make qualitative note of the impact of variation of drop 
size and/or the heat transfer coefficient. The variation 
in time to spontaneous nucleation illustrated in Figure 
5.25, when the drop size is changed, is relatively small. 
This leads to the conclusion that small changes in drop 
size will lead to a minor fluctuation in the time of the 
occurrence of the spontaneous nucleation in the Na, and 
will alter the time scale of the stress plot. However, 
as the induced stresses in the shell are directly related 
to the temperature difference across the solidifying UO,, 
the magnitude of these stresses will not significantly 
change. When considering variation in the heat transfer 
Boeetfticient, this is no longer true, as the interface 
temperature will change. But, as illustrated in 
Section 5.6, when h is varied, the stress levels are still 
high enough to result in fracture. Qualitatively then, 
varying h will alter both the stress level and the time 
to spontaneous nucleation, but not the outcome of the overall 
interaction. 

In conclusion, based on this basic first-order 
analysis, it appears reasonable to assume that the only 
constraint around the Na particle is that provided by the 
molten and partially solidified (but not contiguous) UO.- 


Therefore, the speculation by Fauske (74) that spontaneous 
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nucleation (due to entrainment, wetting and superheating) 
leading to a small-mass vapor explosion is the fragmentation 


mechanism for Na trapped in UO. appears to be a plausible 


2 


description of the interaction. 





CHAPTER SIX 


THERMAL STRESS FRAGMENTATION IN ALTERNATIVE 


MATERIAL COMBINATIONS 


mere tntroduction 
Qualitatively, experimental results of FCI 
simulations using Alumina (A103) as the hot material 
interacting with a sodium coolant (126) have shown 
behavior similar to that observed in UO.-Na interactions. 
The interactions were relatively mild and no strong shock 
waves were observed. Also, the degree of fragmentation of 
the molten Alumina was extremely high. A secondary 
experiment done with Alumina shock-cooled in an Argon 
atmosphere exhibited fragmentation nearly identical to 
that observed in the A1,0,-Na experiments. The 
Similarity of the materials and the resultant interaction 
lead one to believe that the Al,0,-Na Simulation is 
comparable to the actual UO,-Na interaction. Further, 
the similarity between the thermally shocked fragmentation 
residue and that observed in the simulation experiment 
seem to lend credence to the thermal stress mechanism 
of fracture under discussion here. Thus, an examination 
Of both the Al.,03 


Varying interface temperature was undertaken to assess 


-Na interaction and that of Al,0, with 


the applicability of the thermal stress fragmentation model. 
149 
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Furthermore, due to the extensive experimentation 
with metal-water simulants, it is desirable to assess the 
applicability of the thermal stress fragmentation model 
to such a case. Thus, an analysis of a Tin-H,0 
interaction from a thermal stress point of view is 


included, as well. 


6.2 Thermal stress fragmentation of alumina (Al,03) 
Initially, the case of alumina dropped iis Sedum 

comparable to experimental conditions described by 

Lazarrus, et al. (127) was considered. A molten A1,0. 

droplet at its melting point (2040°C) was assumed 

instantaneously immersed in a pool of sodium at 450°C. 

Assuming perfect wetting, the contact interface 

temperature for this case is 1033°C. Based on a fixed 

surface temperature condition, the subsequently developed 

stresses at the external surface of the shell are as 

depicted in Figure 6.la (excluding density change effects). 
The values of Ky developed as the generated crack 

moves across the shell, for the stress conditions at various 

times is illustrated in Figure 6.2. The analysis is 

based on the progressive crack growth model described 

in the last chapter for UO./Na, and it can be seen that 

the K. values are substantially larger than the material 


I 


Koc EOr A103. The value of Koc Ler Al,0., is based on 


the experimental analysis of Davidge and Tappin (127) 


yielding a value for the surface energy at fracture of 
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Figure 6.1b--Tangential surface stress for A103 
with Ty = 1838°C. 





152 


%efE0I = “L 
wo 29'0 

ON - 
O IW 


SNOTIPA 3e YQMOIH YyOeAD satssarzboad Ao;Z 


0'l 





“ON UT fo fiy OF + 


49/9 HLONST WOVYD TIWNOISNSWIGNON 
80 





" * UOIIPUOD edDjuNS 


-- snippy doug 
+ + yupo09 


|OIJOLDW jen 


sdoq4s owt 
M JO eNTeA--7°9 eANHTY 


90 vO 20 


=== = [= a So ee 


SS3NMOIHL TI3HS = 49 


822000 


G2¢~c0'0 


SO8s0'0 


lOocel'O 


(988 P6SI2'O = I 


ae 


00 


Ol 


(NIA ISd 401) ~ |» 





153 


= 6 j/m, and a fracture stress in three point bending 
of 1.595 x 10 psi for an anticipated flaw size of 
100 um. Utilizing the analysis described in Ref. 110, the 
mean value of K was determined to be 1.5 x mo psi/in. 


re 


The generated K. values based on the total generated 


if 
stress would lead to almost immediate propagation through 
the shell, relieving the internal pressure. This fact, 
coupled with the knowledge that alumina decreases in 
volume by approximately 20% upon solidification (128), 
also yielding reduced internal pressures, necessitates 

an evaluation only considering the thermal stresses. The 
values of K. based only on thermal stresses are shown in 


if 
Figure 6.3. Again, it can be seen that Koc is substantially 
exceeded. This fact is more easily seen in the expanded 
view of Figure 6.3. Here, the value of Koo is plotted 
against the actual, and not the non-dimensionalized, 
crack depth. The interesting aspect of this figure is that 
the values of K. for all times are nearly on top of one 
enother. It is of further significance to note that 
Koc is exceeded for flaw sizes less than the anticipated 
100 wm in all cases. Therefore, if the actual flaw size 
is around 100 um, the crack movement can be anticipated 
as being initiated with an excess of stored energy, 
leading to a propagation through the shell. Thus, as in 


the case of UO./Na, the thermal stresses alone appear 


Significant enough to result in complete shell failure. 
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Figure 6.3--Value of K. for progressive crack growth at 
various tithe steps for A150, in Na (thermal 
stress only). 
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Figure 6.4--K_ versus actual crack depth for Al 9) 3 
(thermal stress only). 
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An analysis based on varying the interface 
temperature to simulate experimental conditions where a 
varying h occurs at the interface was carried out for Al,0.. 
The results for an interface temperature of 1838°C are 
milustrated in Figures 6.1b, 6.5-6.7. Again it is seen 
that flaws greater than 100 um will be propagated by 
thermal stresses alone. However, when the interface 
temperature 1S increased to 2020°C the required flaw 
size increases significantly, and greatly exceeds the 
100 um size. This fact is graphically illustrated in 
Figure 6.8. As the surface temperature will decrease with 
time, the flaw size required for fracture will also 
decrease. Thus, as the required flaw size approaches the 
actual flaw size, fracture will occur. In the case of 
thermal stress fracture in such a simulation, the delay | 
time will be related to the changing interface 
temperature rather than the time to develop the shell and 
generate the stresses (the case for Al,0, in Na). The 
fracture mechanism, however, will be the same in both 
cases. 

These calculations appear to verify the experimental 
observations of Ref. 126 and further substantiate the 
thermal stress mode of fragmentation for brittle ceramic 
materials. Furthermore, the experimental simulation of the 


U0,/Na interaction appears to be well modeled through the 


use of Al,0,- 
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Figure 6.6--Value of K, for progressive crack growth 
at various time steps for Al,0, (thermal 
stress only) for T, = 1838°C” 
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6.3 Applicability to the tin/H.0 interaction 

The numerous experimental cd mmprciene Mite dM 229811 6) 
molten metals dropped into water, many involving tin/H.,0 
interactions, led to a desire to determine the applicability 
of the thermal stress fragmentation mechanism to these 
Simulations. It 1S apparent from the outset that molten 
metals heated significantly above their melting point 
would invalidate the solidification phenomenon over the 
time scale of interest. This would immediately invalidate 
the thermal stress fragmentation model for many of the 
existing high temperature experimental results. As well, 
the extreme ductility of metals such as lead, aluminum, 
and tin led to some doubt as to the applicability of 
linear-elastic fracture mechanics and brittle fracture 
mechanisms. In light of these observations, an estimation 
of the plastic zone size expected in the case of tin/H.0 
(for tin at temperatures not substantially in excess of 
the melting point; 232°C) was carried out. 

In the case of tin/H.0O, the interface temperature is 
known to vary with time, as in the case of Al.0,/argon. 
Thus, in an attempt to gain some insight into whether the 
thermal stress fragmentation model was applicable or not, 
the interface temperature was varied from 200°-100°C, 
while employing the fixed surface temperature solution. 
The values of Ky developed in the shell for relatively 


comparable times based on differing interface temperatures 
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are illustrated in Figure 6.9. Even though the volumetric 
reduction upon solidification in tin is only 2.6%, in an 
effort to make the analysis conservative, the pressurization 
was neglected and the resulting Ky values are seen in 

Figure 6.10. 


i@estadius of the plastic zone around the tip of a 


flaw in the shell can be estimated as (129) 


me | Sy 
ees ae 
N4 


where ag 1s the tensile yield stress. For tin, a mean 


yae 
CG7-1)) 





value at elevated temperatures of 800 psi was selected 
fOr the yield stress (132). Figure 6.11 illustrates the 
variation of the size of the induced plastic zone with 
increasing flaw size. It can be seen that in all cases, 
for flaws less than 10% of the shell thickness in size, 
the plastic zone diameter exceeds the shell thickness. 
This percentage rapidly decreases with decreasing 
interface temperature. For such a case, KS loses its 
validity and complete plastic deformation will take place 
Over the tensile zone of the shell (130). Another way of 
looking at this phenomenon is through the gross approximation 
that the behavior will be elastic as long as the plastic 
zone size is on the order of or less than the crack 
depth. This would only apply to small cracks, as the 
Plastic zone size still must remain significantly smaller 


than the specimen dimensions. A line representative of this 
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Figure 6.9--K_ values for progressive crack growth at 
various interface temperatures for tin in 


HO. 
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Figure 6.10--K. values for progressive crack growth 
at various interface temperatures for 


fam on HO (thermal stress only). 
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Figure 6.11--Estimation of generated plastic zone 
size versus crack depth for molten tin 


in H,O at various interface temperatures. 
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criterion is also included in Figure 6.11 for crack 
depths less than 10% of the shell width. Based on this 
analysis, the model of thermal stress fracture would not 
be applicable to ductile materials such as tin, lead, or 
aluminum, and an alternative mechanism must be considered 


responsible for the fragmentation in these materials. 





CHAPTER SEVEN 
CONCLUSION 


ml summary 

In this analysis, it has been determined that during 
a U0,/Na thermal interaction the stresses generated in the 
solidifying UO, are sufficient enough to result in fracturing 
of the solid shell. This stands then as analytical 
Werification of the original postulations of Cronenberg 
et al. (77). The value of Ky resulting from the stress 
field in the shell exceeds the critical fracture 


toughness, K of UO. by such a margin that rapid, 


ne” 
instantaneous crack propagation is anticipated. It has 

been determined that the high initial Ky values will result 
in the opening of an initial fissure, thereby relieving 

the internal pressure. In this case, any subsequent fracture 
would have to be the result of thermal stresses alone. When 
only considering the thermal stresses, the generated value 

of Ky was still found to be sufficient to cause fracture 
initiation. Furthermore, a qualitiative assessment of the 
crack growth stability led to the determination that even 
though Ky is driven to zero within the shell, due to the 


compressive stresses, the dynamic crack growth should be 


enough to allow propagation through the shell. This would 
167 
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maen allow for multiple cracking from other flaws, even 
after the initial one relieves the pressure. 

If, perchance, the initial pressure-relieving fissure 
fails to develop, the density change during solidification 
would cause an internal void to develop within the 
solidifying droplet, reducing the pressure. The actual 
pressure level existing within the shell is in some doubt, 
in this case. It could be as low as the saturation pressure 
of uo. Or, if gases have been driven into solution during 
the initial stages of solidification while the drop is 
pressurized, it could be some combination of the 
non-condensable gas pressure and the Saturation pressure, 
as the non-condensables come out of solution during void 
Formation. In any case, during the initial stages of 
solidification the pressurization is present and the 
stresses are high enough to result in immediate fracture. 

In a separate parametric study, it was found that 
variation of the surface heat transfer coefficient did not 
Significantly affect the conclusions drawn from the fixed 
Surface temperature solution. The surface heat transfer 
coefficient was decreased from De aes - Salil ting 
(comparable to T,-constant) and the total stress level was 
found to increase. This was primarily a result of increased 
pressures working over a thinner shell, as solidification 
was retarded, raising the pressure stress. The thermal 


stress, however, was found to decrease. Thus, initial 
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fracture under pressurization is more likely. Here again, 
even with the relief of the internal pressure, the thermal 
stresses appear significant enough to cause fracture. 

In the inverse problem, when small amounts of liquid 
sodium become entrapped in larger amounts of molten UO., 
it was determined that the total stresses were high enough 
to result in the fracture of the solidifying UO,. Asa 
result, the only constraint estimated as existing by the 
time the trapped sodium has been heated to its spontaneous 
nucleation temperature is the surrounding partially 
solidified and molten UO... Therefore, the spontaneous 
nucleation theory, uninhibited by mechanical constraint, 
appears to be directly applicable in this case. 

Extending the analysis to experimental simulants led 
initially to an investigation of themal stress fracture of 
A150, in sodium. The analysis demonstrated that the 


behavior of A1,0 was the same as that for U0,, with 


5 
fracture imminent under the generated stresses. As well, 
the volume change on PT dirieacion is similar, further 
werafrying comparability. In an effort to ascertain the 
Beasibility of utilizing Al1,0, for experimental verification 
of the thermal stress fragmentation model, the surface 
temperature was varied to simulate a range considered 
applicable to a thermal quench. The results indicate 


that initially, while the surface temperature is still at 


a high level, solidification will proceed but fracture 
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will not initiate, as the critical flaw size required is 
‘excesSive. As the surface cools, however, Ky will 
increase (decreasing the critical flaw size) and fracture 
will follow. Even though the perfect wetting and fixed 
surface temperature conditions which are characteristic 
in sodium may not be present in this case, the fracture 
mechanism and fracture behavior of the A140, would be 
Similar to that occurring in a sodium quench. The major 
difference would be the delay time to fracture. 

An evaluation of the potential applicability of the 
thermal stress fragmentation model to metal/water interactions 
was conducted as well. In this case, a Significant factor 
came to light. That is, that in the case of ductile 
materials such as lead, tin, and aluminum, the plastic zone 
Size was found to grow to sizes approaching the specimen 
dimensions for very small flaws, at short times. This 
results in invalidating the critical stress intensity 
factor approach, as gross plastic yielding would take 


place in the tensile zone of the solidifying shell. This 


fact, coupled with the knowledge that metal/water 


interactions result in fragmentation even when the molten 


material temperature is high enough to prevent solidification 


over the time span of interest, has led to the conclusion 


that fragmentation in such a system is the result of an 


alternative mechanism. The validity of the thermal stress 


fragmentation model is, therefore, limited to semi-brittle 
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BreOrittle materials which, for the most part, behave 
elastically during the interaction process. 

The overall conclusion reached is that thermal stress 
initiated fracture is a viable description of the 
fragmentation process in the UO./Na System. Furthermore, 
it 1S concluded that other, easier to handle brittle 
ceramics, such as Al1.,0,, will sufficiently model the 


interaction to allow correlation of theory and experiment. 


7.2 Recommendations for future work 
7.2.1 Experimental 

From the standpoint of experimental study, it is highly 
recommended that a thermal shock-type experiment 
mei lizing Al,0, be carried out. This would allow 
verification and correlation of experimental and theoretical 
work. An example of the types of correlations which would 
be meaningful is illustrated by Figures 7.1 and 7.2. These 
figures are based on the limited data generated in this 
Study and are not all-inclusive. Figure 7.1 illustrates the 
variation of Ky with time for various values of the 
interface temperature, based on the most probable flaw size 
of 100 um. Experimental values would be directly obtainable 
for comparison on such a basis. The interesting aspect of 
this figure is that it demonstrates the fact that fora 
Specific flaw size, an interface temperature cutoff exists 
above which fracture would not occur. Eventually, 


determination and verification of the value of this cutoff 
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Figure 7.2--K. versus flaw size at various interface 
témperatures. 
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temperature for the UO0,/Na system should be carried out, 

as it would seem to indicate that if the coolant 

temperature reached a certain level all further fragmentation 
would be suppressed. Figure 7.2 provides for a slightly 
different point of view which is equally applicable to 
experimental evaluation. In this case, if the time span 
envelopes are expanded analytically to cover the times 


rom t. 
m 


in (governed by solidification kinetics) to t. 


ax 
(complete solidification), the anticipated time to fracture 
for various interface temperatures could be predicted and 
subsequently verified experimentally. Note that the cutoff 
temperature phenomenon is exhibited in this figure for a 
100 um flaw size, as well. 

In evaluating coolant materials for experimental use, 
it was originally considered desirable to use an inert gas 
Such as argon, as it would help to eliminate complications 
from second-order effects (such as boiling, in the case of 
water) while still allowing for "observability". However, 
in conducting preliminary heat transfer analyses, it was 
found that the poor heat transfer characteristics of gases 
required either a tremendously long fall, for dropping 
experiments, or an extremely high velocity for injection 
experiments to allow the interface temperature to reach 
levels capable of generating significant stress levels. 
Investigation of other feasible materials led to the 


consideration of molten salt (NaCl). It was found that 
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for a molten salt coolant at approximately 900°C, the 
interface temperature would be approximately 1074°C. This 
value of the interface temperature is sufficient to result 
in fracture within reasonable times as compared to gases. 
The adequacy of interface temperatures in this range was 


demonstrated in the earlier analysis of hte Al 0, interaction 


2 
@nd in Figures 7.1 and 7.2. Thus, utilization of a molten 
salt coolant would appear to provide a reasonable experimental 
evaluation of the interaction. An experimental verification 


of this type is considered essential to the final acceptance 


of the thermal stress fragmentation mechanism. 


Meee2 Analytic confirmation 

In parallel with the above experimental work, a 
verification of the results obtained thus far for Al,03, 
based on a more streamlined calculation procedure, is 
desirable. Specifically, a numerical solution, such as 
that of Tao (134) would be the most useful, as it provides 
the capability of including various heat transfer 
coefficients at the surface, if required. Subsequently, 
further improvements such as including the finite element 
solution for predicting K, values would be beneficial. 
Then, once an accurate prediction of the surface heat 
transfer condition/coefficient is obtained, theoretical 
predictions of increased accuracy and reliability would be 


available for comparative purposes. 
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It would also be desirable to make a more detailed 
analytic analysis of the crack growth and stability under 
the influence of thermal stresses alone, thus making it 
possible to analytically verify what are at present 
gualitative observations. Further, the anticipated 
inherent flaw size and distribution require more accurate 
definition. This would then allow determination of the 
number of initiating cracks, and the resultant particle 
Size distribution, with the possibility of further 
extension to include variation with time. 

One other aspect, of higher order, which requires 
analytical interpretation concerns the actual, 
nonsymmetric solidification process. The density change, 
Besulting in void formation, will lead to a region in 
which there are significant temperature and stress field 
perturbations. The impact of these effects, anticipated 
as being small, awaits detailed analysis. 

These analytical improvements to the thermal stress 
fragmentation model, as it stands, coupled with the 
experimental work outlined above would surely provide a 
strong foundation for the concept, as a whole. The work 
accomplished to date is only the first step toward this 


end. 
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APPENDIX SC 
ANALYTIC SOLUTION TECHNIQUE 


The purpose of this appendix is to summarize and 
outline the basic solution technique employed when using 
Meechnermal stress fracture model. This, hopefully, will 
allow future workers to utilize or to extend this model 
without having to start from the beginning. 

The primary tool employed is the computer solution, 
based on the fixed surface temperature condition, included 
in this appendix. The thermal solution employed is that 
of Adams (82) which, through numerical integration, 

Me vides the times for successive positions of the 
solidification front, and the temperature profile in the 
shell at these times. Subsequently, this information is 
used as input to equations 3.1-3.3 which are also 
Wmerically integrated to yield the induced pressure and 
stresses within the shell. Note that the pressure here is 
a result of a balance between the attempted thermal 
Meatraction of the shell and the slightly compressible 


‘molten UO The primary inputs, variable designation 


2° 
and the output of the program are all included on comment 
cards in the program listing included herein. Operation 


of the program is self-explanatory. Figures 5.10-5.12 
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are illustrative of the type of information obtained 

from the program. As well, the complete shell stress 
distribution for a given time is obtained, as seen in 
Mmogures 5.18 and 5.19. 

For convenience in checking the program, a sample 
O@meput for 2800°C uO, in 200° sodium with the interface 
temperature at 945°C is included following the program 
listing in this appendix. This was a run for a 
solidification front position of .3467 cm and it employed 
the vo. property values given in Appendix A. (Note: 
ensure that the thermal contraction coefficient value used 
is converted to "per °C" to be compatible with the rest 
of the input data.) 

Utilization of the stress distribution at various 
times (e.g., Figure 5.19) and equation 5.17 allow 
determination of K, for various crack depths as illustrated 
in Figures 5.20 (for total stress) and 5.21 (for thermal 
Stress only). A comparison with the material fracture 


toughness, K then allows determination of whether 


ine 
fracture will occur or not. Finally, an investigation of 


the point where K, exceeds Koc for a known flaw size leads 


> 


to determination of an anticipated delay time to fracture 
for the material under consideration. 

Meee that a critical input to this solution is the 
surface temperature of the molten drop, and a determination 


of its value must be made prior to initiating any analysis. 





eZ 


There iS a provision for varying the surface temperature 

in the program to allow for parametric evaluations which 
simulate the surface temperature variation anticipated 

when a surface heat transfer coefficient exists. This 1s 
only good for cursory estimates as were made in this 
report, as the time rate of variation is a critical factor 
in any meaningful analysis. Thus, a solution such as that 
of Tao (134) is recommended when it is known that the fixed 


surface temperature condition is not applicable. 
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